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This set of four lectures reviews various aspects of the theory of a dilute low-temperature trapped 
ON ' Bose gas, starting with (I) a review of the Bogoliubov description of the elementary excitations 

ON ' in a uniform system. The treatment is then generalized (II) to include the new physical effects 

of a confining harmonic trap potential on the condensate and its normal modes. An equivalent 
^ ' hydrodynamic description (III) focuses directly on the density and velocity fluctuations. The physics 

of vortices (IV) in an incompressible fluid is summarized and extended to the case of a trapped Bose 
■ condensate. 
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I. UNIFORM DILUTE BOSE GAS 



Bogoliubov's Q seminal 1947 paper originated the modern theoretical description of a dilute interacting low- 
temperature Bose gas. He considered how a weak repulsive interparticle potential affects the ground state and low- 
lying excited states of a uniform condensed Bose gas, showing that the interactions qualitatively alter the dispersion 
"olj ' relation at long wavelengths from the familiar quadratic free-particle form p 2 /2M to a linear "phonon"-like structure 
sp with speed of sound s. 
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A. Brief review of a uniform ideal Bose gas 

The standard textbook example |2|| of an ideal Bose gas treats a uniform system of N particles in a cubical box 
of volume V — L 3 with periodic boundary conditions and number density n — N/V. The relevant single-particle 
l— "™ '. states are plane waves F -1 / 2 exp(ik • r), with k = (2n/L)(ni, n 2 , n^), where rii is any integer; the energy of such a 
' state is e" = Ti 2 k 2 /2M, where M is the particle's mass. In the classical limit, the thermal de Broglie wavelength 
J> • At = (2irh /MhsT) 1 / 2 is much smaller than the interparticle spacing I w n -1 / 3 , so that diffraction effects are 
' negligible. The inequality At -C I necessarily fails as the temperature is reduced (or the number density increases), 
and quantum diffraction becomes important at a temperature T c determined by the approximate condition At c ~ I 
' (or, equivalently, nA Tc w 1). A detailed calculation yields the exact expression 

5o ■ nA l = C(l) * 2 -612, (1.1) 

"j***; i where C(|) is a Riemann zeta function. The transition at T c represents the onset of the special form of quantum degen- 
eracy known as Bose-Einstein condensation. Alternatively, the transition temperature T c is given by the approximate 
£h , relation 

13 1 h 2 n 2 / 3 
2 ' k B T c a ^— (1.2) 

g . M v ' 

Q ' that characterizes the onset of degeneracy in any ideal quantum gas, providing a qualitative description of electrons 
in metals and white dwarf stars, nucleons in nuclear matter and neutron stars, and the fermionic isotope 3 He, as well 
as the bosonic isotope 4 He. 

For T < T c , a uniform ideal Bose gas has a macroscopic number of particles Nq(T) in the lowest single-particle 
jj^ ■ state (here, that with k = 0), given by 

_ i 

Nq(T) / T \3/2 

- 1 - — for T < T c . (1.3) 



N \T C , 

Thus a finite fraction of all the particles in a condensed Bose gas occupies a single quantum-mechanical state. In the 
limit T = 0, this fraction is unity; all the particles in a uniform ideal Bose gas at zero temperature and fixed volume 
V have zero momentum (and therefore exert no pressure on the confining walls). 
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B. Effect of weak repulsive interactions 

Bogoliubov [H introduced the microscopic treatment of a "weakly interacting" or "nearly ideal" Bose gas, starting 
from the second-quantized hamiltonian 

H = e fc a ic a k + ^7 Vk!-k3 a ki a k2 a k3 a k4 <5 kl+k2ik3+k4 , (1.4) 

k ki,k2,k3,k4 

where the Kroneckcr 5 ensures momentum conservation, and — J d 3 r exp(— zk • r)V(r) is the Fourier transform of 
the interparticle potential. The operators a k and a k obey bosonic commutation relations 

[a k ,a k ,] = 5 kik /, [a k ,a k '] = [aL a k'] = °- ( L5 ) 

Bogoliubov's basic idea is very simple and elegant. If the interacting Bose gas is sufficiently dilute, it should differ 
only slightly from an ideal gas, so that most of the particles in the true interacting ground state should have zero 
momentum. Furthermore, only two-body collisions with small momentum transfer are significant, and these may 
be characterized by a single parameter known as the s-wave scattering length a (here generally taken as positive, 
corresponding to a repulsive interaction). For a strong repulsive potential, a is simply the range, and the more general 
case is treated below. In the dilute limit a <C I (or, equivalently, na 3 <C 1), the interparticle potential may be 
approximated by a "pseudopotential" V(r) w gS^(r), with a constant Fourier transform Vk = g, and the original 
hamiltonian in eq. ( [l.4| ) becomes 

H = X/° a k a k + 2V al i al 2 a^ 3 a^S iil+i , 2M+k4 . (1.6) 

k ki,k2,k 3 ,k4 



C. Review of scattering theory 

The model interaction strength g can be related to the more physical parameter a by requiring that the many- 
body hamiltonian (1.6) reproduce the correct two-body scattering in vacuum fl. The Schrodinger equation for the 
scattering wave function V'k(r) of two particles with mass M and initial relative wave vector k can be rewritten in 
terms of the relative separation r = ri — V2 and the reduced mass \M 

- |^ k (r) + ^(r)Vk(r) = ^k(r), (1.7) 

where E = 2e° = H 2 k 2 /M. With the outgoing-wave Green's function of the Helmholtz equation [Q, this partial 
differential equation can be recast as an integral equation 

yr f p ih\r— r'| 

^ k(r) = e ' k ' r ~ J ^ "FT rT y(r ' } Mr ' 1 (L8) 



For simplicity, assume that the potential has a finite range (and no bound states), so that eq. (1.8) has the asymptotic 
form for Irl = r — * oo 



e lkr 



<Mr)~e 4k - r + /(k',k)— . (1.9) 

This expression is the sum of an incident plane wave exp(ik • r) with wave vector k, and an outgoing scattered spherical 
wave r" 1 exp(ifcr) with a proportionality factor 

/(k', k) = - J|y j d\' e- lk ' r ' V(r>) Vk(r') (1.10) 

known as the scattering amplitude for a transition from the initial relative wave vector k to a final relative wave vector 
k'. Note that this definition remains valid even for a singular potential because it involves the true wave function tp-^ 
that vanishes wherever V(r) diverges. In contrast, the Born approximation takes ^k(r) ~ V'k ( r ) = ex p(*k • r), and 
the resulting Born approximation to the scattering amplitude 
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/ B (k',k) = --^y k ,_ k (1.11) 

Ann 

clearly requires that V(r) be integrable and hence non-singular. 

For a spherically symmetric potential with k' = k (so that energy is conserved), the scattering amplitude depends 
on k and the scattering angle 9 (defined by cos 8 = k ■ k ) 

/(M) = X)— r^e* sin<S,P z (cos0), (1.12) 

where Si is the energy-dependent phase shift for the Zth partial wave. In the low-energy limit (k — > 0), the s-wave 
phase shift reduces to <5o ~ —ka, which defines the s-wave scattering length a, and higher-order phase shifts are of 
order k 2l+1 for I > 1. Thus, only the s-wave term (I = 0) contributes as k — > 0, and the scattering amplitude takes 
the simple form 

lim /(k',k)«-a. (1.13) 

|k| = |k'H0 

For many purposes, it is convenient to introduce a Fourier decomposition of the scattering wave function, with 

M') = J (§? eip ' r ^(p), (1-14) 



and the exact scattering amplitude in eq. (1.1C) can be rewritten as a convolution of the Fourier transforms of the 
interaction potential and the wave function 



/(k ' k) ^-^/|^ fk '- pV:k(p) - (L15) 



In addition, the outgoing-wave Green's function has the Fourier representation 

e *|r-r'| f d 3 giq-(r-r') 



r — r'| J (2ir) 3 k 2 — q 2 + i-q 1 



(1.16) 



where r/ — > + defines the appropriate contour in the complex q plane. A combination with the integral equation (l.S) 
then yields an exact expression for the Fourier transform "0k (p) 

M P ) = (2nf ^)(k - p) - j"_ /( P' k L (1.17) 



incident wave 



k 2 — p 2 + IT] 



cattcred wave 



again expressed as the sum of an incident wav e and an ou tgoing scattered wave. 

As a final step, a combination of eqs. ( 1.15| ) and (1.17) gives a corresponding integral equation for the scattering 
amplitude 

4 * B W) = v,_ k -4^^%^i>, (L18) 



M J (2tt) 3 k 2 -p 2 +i v ' 

which requires kno wledg e of the scattering amplitude for p 2 ^ k 2 (known as "off-the-energy-shell" ) . Note that the 
leading term of eq. (1.18) is simply the Born approximation fg, and an iteration gives the first correction 

/(k ,k) « y k ,_ k + ^ / —-3 (1-19) 



M h 2 J (2tt) 3 P - p 2 + ir? 

In the limit fc' 2 = fc 2 — > 0, this equation becomes a perturbation expansion for the scattering length 

^--^/^ + (L20) 



M ?i 2 7 (2tt) 3 p 2 

and the first correction is finite if V p vanishes sufficiently rapidly as |p 
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The situation is different in the singular case of the pseudopotential V p = g, for eq. ( 1.20| ) reduces to 



4nh 2 a Mg 2 f d 3 p 1 , 

~9--t 7^3 + ---> (1-21) 



(1.22) 



M " h 2 J (2tt)3 p 
To leading order, this equation relates the pseudopotential g to the measurable scattering amplitude 

Ann 2 a 



It typifies the general prescription that a first-order perturbation calculation can usually be extended to treat a dilute 
hard-core gas by substituting the true s-wave scattering length a for the first Born approximation ag — MV / 4nh 2 ; 
in the present context of a dilute Bose gas, thi s pro cedure was suggested by Landau (see last footnote of ref. 0). 



In contrast, the second-order correction in eq. ( |l.2l[ ) diverges for large values of |p|. As seen below, the Bogoliubov 
theory yields a similar divergence for the ground-state energy that disappears when re-expressed in terms of the 
physical quantity a (for an alternative treatment of the divergence arising from the use of the pseudopotential, see 
ref. |). _ 

There is one remaining subtlety, for the s-wave scattering length has been defined as if the particles were distin- 
guishable, whereas the true overall wave function is symmetric. Consequently, the actual differential cross section is 
obtained from a symmetrized scattering amplitude 

^ = \.f(k,9) + f(k,7T-e)\ 2 , (1.23) 

which reduces to 1 2a. | 2 = 4a 2 in the low-energy limit, four times that for distinguishable particles. The corresponding 
total cross section o~t = 87ra 2 is obtained by integrating over one hemisphere, because the particles are indistinguish- 
able. 



D. Bogoliubov quasiparticles 



Given the approximate relation in eq. ( 1.22 ) between the pseudopotential g and the s- wav e scattering length a, it is 
now feasible to proceed with Bogoliubov s jl] analysis of the model hamiltonian in eq. (1.6), following the treatment 



of ref. |J, sees. 35 and 55. The basic observation is that an imperfect Bose gas at low temperature should resemble 
an ideal Bose gas to the extent that a macroscopic number Nq of particles occupies one single-particle state. In 
the simplest case of a stationary uniform system, this condensation occurs in the mode with k = 0, but many 
other possibilities also can occur (for example, a condensate moving with velocity v). The effect of the repulsive 
interactions is to scatter particles from the ideal-gas ground state k = to pairs of higher-momentum states with 
±k 7^ 0. This depletion of the zero-momentum condensate implies that Nq is definitely less than N, even at T = 0, 
but the "condensate fraction" Nq/N is assumed to remain finite in the thermodynamic limit (N — > oo, V — » oo, with 
N/V fixed); this assumption must naturally be verified at the end of the calculation. 

Recall that the ground state of an ideal uniform Bose gas \$o(N)) = \N, 0, 0, • • •) = (aJ) Ar (iV!)" :L / 2 |0, 0, 0, • • •) has 
all iV particles in the single-particle mode with k — 0. The associated annihilation and creation operators then yield 
macroscopic coefficients when applied to this ground state 

a \$ (N)) = VN\$ (N-l)), al\$ (N)) = Vim \$ Q (N + 1)) , (1.24) 
whereas their commutator 

[a o ,<4]|$o(iV)> = |$o(A0) (1-25) 

gives a coefficient unity. To the extent that the thermodynamic properties of an ideal Bose gas are independent of 
the addition or subtraction of one particle, it is natural to treat the condensed mode differently from all the others, 
replacing the operators ao and a by pure numbers, because their commutator is of order N~ 1 / 2 relative to their 
individual matrix elements. 

The same situation is assumed to describe the ground state |$) of an imperfect uniform Bose gas, as long as the 
condensate density 

<Mpm = ^ sna (L26) 
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remains finite in the thermodynamic limit. Consequently, Bogoliubov proposed the now famous prescription 



aa, a 



'N , 



(1.27) 



with the operators ao and aj reinterpreted as pure numbers that commute. For simplicity, he also imposed the 
additional stronger condition of small total depletion N — No <C N, which limits the strength of the interparticle 
potential and requires that the temperature T be small compared to the critical temperature T c for the onset of Bose 
con densa tion. It is important to realize that this second assumption is distinct from the Bogoliubov prescription in 
eq. (1.27), which merely requires the existence of a finite condensate density no in the thermodynamic limit (the more 
general case at finite temperature is descri bed, f or example, in refs. as we ll as in Griffin's lectures in this volume). 
Substitute the Bogoliubov prescription (1.27) into the model hamiltonian (1.6) and keep only the leading terms of 



order Nq and A^o (for a uniform gas, the terms of order N vanish identically because of momentum conservation 
but a similar situation holds in general). The resulting interaction part of the hamiltonian becomes 



TT ^ 

int ~ 2V 



Nl + 2N ( a W + a-k«-k) + N oJ2 ( a k a -k + ft ka-k 



k^O 



k^O 



(1.28) 



In the present context, it is convenient to consider only states with a fixed number N of particles (see sec. [n] for a 
more general approach using the grand canonical ensemble). The number operator 



N — a ic a k ~ ^o + \ E ( a ic a k + aL k a-k) 



(1.29) 



k^O 



can then be replaced by its eigenvalue N, and eq. ( 1.29|) t hus serves to eliminate No in favor of N through terms of 
order iV 2 and N. As a result, the model hamiltonian (1.6) becomes 



H 



gJP 
2V 



ng) \ala k + a^ k a- k ) + ng ( a k a_ k + a k a_ k 



k^O 



(1.30) 



This approximation evidently neglects terms like (X^k^o a iL a k) 2 j which are small for N — No *C AT. The Bogoliubov 
hamiltonian (1.30) is a quadratic form in the creation and annihilation operators and can be diagonalized with a 
canonical transformation. 

To demonstrate this remarkable feature, Bogoliubov introduced a new set of "quasiparticle" operators a k and a k , 
defined by the linear transformation 



a k = u fc a k - v k a f _ k 



(1.31a) 



,t _ 



u <fc a -k - v k a k , 



(1.31b) 



with real isotropic coefficients u k and v k . Note that both operators a k and a_ k reduce the total momentum by fik. 
This transformation is canonical if the new operators also obey bosonic commutation relations, with ^,0^] = 5k, k' 
and [ak, a k <] = [a k , a k ,] = 0. It is easy to verify that this condition requires 



u\ - v\ = 1 for all k ^ 0. 
Direct substitution of eq. ( 1.3l|) into the Bogoliubov hamiltonian ( 1.3C ) yields 

H = \gn 2 V + ^2 [( e k + n 9) v k - ngu k v k ] 



(1.32) 



k^O 



k^O 



I E \- ng + v k) ~ 2 ( e fe + n 3) u k v k] (a k a -k + «ka_k 



(1.33) 



k^O 
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Here, the first line is a pure number with no operator character, and the second line is diagonal in the quasiparticle 
number operator a^a^. In contrast, the third is off-diagonal in the quasiparticle number, but it can be eliminated 
entirely by choosing the coefficients Uk and Vk to satisfy 



ng (u 2 k + v\) = 2 (eg + ng) u k v k . 



The subsidiary condition u 2 , — v\ = 1 in eq. (1.32) can be incorporated automatically by writing 

Uk = cosh 9 k and Vk = sinh 9k , 



and a combination with eq. (1.34) readily yields 

tanh 29 k 



ng 



(1.34) 
(1.35) 
(1.36) 



If g is negative, the right-hand side of this equation will diverge for some k, implying unphysical values for a range of 
the 9^ s. The present case of a uniform dilute Bose condensate therefore requires repulsive interactions, with both g 
and a positive. 



The resulting equations for u\ and vi are easily solved to yield 



2 2 

v k = u k 



ng 



E k 



(1.37) 



where 



E k = ^{el+ngf -{ngf = ^( e °) 2 + 2nge% (1.38) 
is defined as the positive square root. Substitution into eq. ( |l.33|) gives the final and simple quasiparticle hamiltonian 



(1.39) 



k#0 



which involves only the quasiparticle number operator for each normal mode k. Since a^a^ has the eigenvalues 
0, 1,2, it follows immediately that the ground state |$) of this model hamiltonian is simply the quasiparticle 
vacuum, defined by 



a k |$)=0 for all k t^O; 



(1.40) 



evidently, |<fr) is a very complicated combination of unperturbed particle states, for neither nor annihilates it. 
In addition, the ground-state energy is given by 



E„ 



\gn 2 



k/O k/O 



(1.41) 



The physics of the quasiparticle hamiltonian is very transparent, for all the excited states correspond to various 
numbers of non-interacting bosonic quasiparticles, each with a wave vector k and an excitation energy 



Ek = 



' ngh 2 k 2 



AI 



4irh z 



h 2 k 2 

2M 



A'kTi an 
M 2 



AI 



hk, for k — > 0, 



for k — > oo. 



(1.42) 



(1.43) 



In the long- wavelength limit (k — > 0), the elementary excitations represent sound waves with the propagation speed 

i ng h 



M M 



V4j, 



(1.44) 



G 



which again shows that g and a must both be positive. For short wavelengths, in contrast, the spectrum is that 
of a free particle, but shifted upward by a constant "optical potential" gn = Airh 2 an/M arising from the Hartree 
interaction with the remaining particles (equivalently, this shift reflects the forward scattering from the background 
medium, producing an effective index of refraction). 

The cross-over between the linear and quadratic regimes occurs when e k /ng 1. This ratio can be rewritten as 
k 2 /%im,a ps 1, suggesting the introduction of a characteristic length 

£=7^= (1-45) 

that is frequently known as the "correlation length" || or the "coherence length;" as shown in sec. ES, however, 
the term "healing length" is usually preferable in the context of a non-uniform dilute trapped Bose gas. Note that 
£ — » oo for a non-interacting gas (a — > 0); in this limit, E k reduces to the quadratic fre e-par ticle spectrum for all 



k, and the linear phonon-like region then vanishes entirely. The energy spectrum in eq. ( 1.42 ) can be rewritten as 
E k /ng = E k M/(4irh 2 an) = fcg J 2 + fc 2 £ 2 , shown in fig. 1 as a function of the dimensionless combination 



The Bogoliubov transformation (1.31) can be inverted to express the quasiparticle operators in terms of the original 
particle operators 

(*k = u k a k + WfeaL k , (1.46a) 



aL k = u k a\ + v k a k . (1.46b) 

The quasiparticle operators are coherent linear superpositions of particle and hole operators, with u k and v k as the 
weight factors. In the phonon regime (fc£ <C 1), these coefficients are both large, with u 2 ~ v\ « (\/8fc£) -1 3> 1, and 
the long-wavelength quasiparticle operators therefore represent a coherent nearly equal admixture of a particle and a 
hole. In the free-particle regime 3> 1), however, the behavior changes to u 2 w 1 and w 2 , s=s (4fe 4 £ 4 ) _1 -C 1, so that 
the short-wavelength quasiparticle creation operator is effectively a pure particle creation operator. Figure 2 shows 
the quantities u 2 and v 2 , as a function of the dimensionless wavenumber k£, illustrating clearly the qualitative effect 
of the repulsive interactions on the long-wavelength properties of a uniform dilute Bose gas. 



E. Macroscopic properties 



Within the Bo goliub ov approximation, the low-temperature behavior can be described with the quasiparticle hamil- 
tonian H in eq. ( 1.39 ) and an associated Gibbs ensemble density operator 



exp(-/? J ff) 



Tr exp(-(3H) 



(1.47) 



where (3 1 = fcgT and the denominator is simply the partition function (the notation Tr denotes a sum over all 
states). Manipulations with this density operator are straightforward, for H is simply a sum of non- interacting 
bosonic contributions, yielding a thermal average over familiar harmonic-oscillator states. For example, consider the 
total particle number at low temperature 



N = N + 



IE 

k#0 



Tr 



(4 



p I a k a k 



,a-k 



(1.48) 



The trace is invariant under a canonical transformation, and it is most readily evaluated in the quasiparticle basis where 



H is diagonal. Use of the Bogoliubov transformations from eq. (1.31) leads to various products of two quasiparticle and 



quasihole operators. Those involving off-diagonal quantities like Tr a k o;l k ^ and Tr (p a k o;_ k ) vanish identically, 
and the diagonal quantities are just the finite-temperature harmonic-oscillator occupation numbers 



Tr(pa k a k ) =f(E k ) 



1 



eM0E k ) - 1 



and Tr (/Sa k a k ) = 1 + f(E k ), 



(1.49) 



involving the familiar Planck (or Bose-Einstein) distribution function f(E k ) for a quasiparticle with energy E k in 
thermal equilibrium at temperature T [note that f(E k ) has no chemical potential because quasiparticles are not 
conserved] . 
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A combination of these results with eq. ( 1 .45 



) immediately yields the relation 



N = N + 



9 



exp(/3£ fe ) - 1 



Nn 



[vl + (u 2 k +vl)f(E k )] 



(1.50) 



Since N is fixed, this equation determines implicitly the temperature-dependent condensate number N (T). In a 
similar way, the low-temperature internal energy E(T) is simply the finite-temperature average of the Bogoliubov 
quasiparticle hamiltonian in eq. (1.39) 



E(T) = Tr (pH) =E g + J2 



k^O 



1 



= E n 



k^O 



(1.51) 



expressed as the zero-temperature ground-state energy E g plus the sum over all thermally weighted excited states. 

For simplicity, consider the total number of non-condensate particles N' = N — N (the depletion of the condensate 
arising from the repulsive interactions). At zero temperature, the thermal factors vanish, yielding 



N'(T = Q)=N- Nq(T = 0) = ^2 

k/O 



(1.52) 



This equation shows that v\ can be interpreted as the number 7V£ of non-condensate particles with wave number k 
in the Bogoliubov ground state |$); since v\ oc (fc£) _1 i s large for fc£ <C 1, the low-lying phonon-like modes have large 
occupation numbers (this behavior reflects the Bose condensation in the mode with k = 0). Nevertheless, the total 
noncondensate number remains finite because v 2 vanishes rapidly for fc£ 3> 1. 



The sum over the excited plane-wave states k in eq. (1.52) can be approximated as an integral X^k^o 
V(2ir)~ 3 J d 3 k ■ ■ •, and the total zero-temperature fractional depletion becomes fjJMfl 



N'(T = 0) N — N (T = 0) 



N 



N 

2na 3 \ 1/2 
n J JO 
8 fna 3 \ 1/2 

3 V 71" 



V dy 



y 2 + i 

(y 4 + ^ 



a/2 



(1.53) 



This result exhibits 



as the relevant expansion parameter, and the assumption of small total depletion requires 



that Vna 3 -C 1. Note that the fractional depletion is non-analytic in the interaction strength a (or g), for the | 
power requires a branch cut in the complex a plane. As a result, this and other corrections to the physical properties 
of an ideal Bose gas are non-perturbative; they cannot be obtained with conventional perturbation theory. Such a 
conclusion here presents no conceptual difficulty, however, for it relies a (non-perturbative) canonical transformation. 
It is not difficult to obtain the first low-temperature contribution to the depletion 



JV'(T) N - Nq(T) 



N 



N 



1/2 



M{k B Tf 
I2h 3 sn 



(1.54) 



It is also interesting to evaluate the ground-state energy from eq. (1.41), which has the form 



E g = \gn 2 V - \ J2 (4 + ng - E k 

k^O 



(1.55) 



An expansion of E k for large fc£ shows that the sum diverges like g 2 X^k^o ^ ' which reflects a failure of second-order 
perturbation theory for a dilute Bose gas. Here, it arises from the assumption that the Fourier transform g of the 
interparticle potential remains constant f or lar ge wave numbers, and the same divergence appeared in eq. (1.21 ) for the 
scattering length to order g 2 . Equation (1.55) can be rewritten by adding and subtracting the divergent second-order 
contribution 



E„ 



hgn 2 v 



12 2 \ ^ 



AI 



k 2 



k^O 



Mg 2 n 2 
h 2 k 2 



E k 



ng 



(1.56) 
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and it is easy to verify tha t the last sum converges. In addition, the first two terms on the right-hand side are precisely 



V times those in eq. ( 1.21 ) for Airh a/M, where a is the physical scattering length. A detailed evaluation yields 



the ground-state energy per particle (|J|] 



E, 



T7= 2 n 5 



Mg 2 



d A k 1 



M 



2i:h an 



(2tt) 3 k 2 
2na^ 1/2 



l2nV 



7T 



v d y 



ng ng 

1/2 



ng 



(/ + 2y 2 )^-y 2 -l + 



1 

V 



M 



128 
15~ 



na 

7T 



1/2 



(1.57) 



which again exhib its the same non-analytic dependence on V na 3 . 

Equation (1.57) expresses the ground-state energy as an explicit function of the number N and volume V, and stan- 
dard thermodynamics provides many other important ground-state properties. For example, the chemical potential 
fj, and pressure p follow from 



and 



P = 



fdEA 
\dNj v 



dE, 
dV 



Anh an 
M 



2ith 2 an 2 



32 (na 



3\V2 



JY 



M 



, 64 /na- 

1 + T {- 



(1.58) 



(1.59) 



In addition, the speed of sound s follows from the compressibility through the general formula s 2 = M 1 (dp/dn), 
giving 



— -V^iran 
M 



3\l/2 



(1.60) 



It is striking that the speed of sound obtained here w ith m acroscopic thermodynamics agrees to leading order with 
the slope of the quasiparticle spectrum given in eq. (1.44); although this zero-temperature relation is not entirely 
obvious, it has, in fact, been verified to all orders in perturbation theory ]9|,|l0|| that the slope of the quasiparticle 
spectrum agrees precisely with the thermodynamic speed of sound. 



F. Moving condensate 



It is not difficult to generalize this analysis to the case of Bose condensation in a mode with non-zero momentum 
?iq, which implies that the condensate moves with velocity v = hq/M. The quasiparticle energy spectrum E^ for 
the state with the proper positive normalization |u| 2 — \v\ 2 — 1 is simply related to the energy E^ for a stationary 
condensate: 

Ek = hk-v + E%. (1.61) 

This expression can be considered a Doppler shift in the frequency; its form agrees precisely with Landau's general 
result in his celebrated 1941 paper on superfluid 4 He |T^,|l2[. In the long-wavelength limit, the excitation energy 
reduces to E^ « hk(v cos 9 + s), where s is the Bogoliubov speed of sound and 9 is the angle between k and the 
flow velocity v. For v < s, the quasiparticle energy is positive for all physical angles 9; if v exceeds s, however, the 
quasiparticle energy becomes neg ative in a cone around the backward direction where cos 9 « — 1. Consequently, the 



quasiparticle hamiltonian in eq. (1.39) becomes unstable, and the system can lower its energy arbitrarily by creating 
quasiparticles with negative energy. This behavior corresponds to the well-known Landau critical velocity v c for the 
onset of dissipation in a superfluid [fni[L2|. In the Bogoliubov model, v c is equal to the speed of sound s. Physically, 



the finite s reflects the presence of repulsive interactions [see eq. (1.44)] and therefore s vanishes for an ideal Bose 



gas. As a corollary, an ideal Bose gas also has v c = and hence cannot sustain superfluid flow. 
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At T = 0, it is straightforward to prove that the total momentum P carried by both the condensed and non- 
condensed particles is P = NMv, where N is the total number of particles, not the condensate number Nq. This 
result implies that the Bose condensation induces a coherence or rigidity, ensuring that, in effect, all the particles 
participate in the motion. The same conclusion also follows with a Galilean transformation from the rest frame of the 
condensate to a frame moving with velocity —v. 

At finite temperature T =/= 0, the corresponding total momentum decreases, according to the expression P = 
(p B /p) NMv, where the reduction factor p s (T)/p < 1 is the temperature-dependent superfluid fraction first calculated 
by Landau |]ll| , |l2| . In his model for superfluidity, the total density p consists of two separate parts, the superfluid 
density p s and the normal fluid density p n , with p s + p n = p. The normal fluid density vanishes at T = and depends 
solely on the excitation spectrum in the rest frame of the condensate. For the Bogoliubov energy Ek , the normal fluid 
fraction p n /p is proportional to T 4 at low temperatures, similar to the phonon contribution in superfluid 4 He [pSf. 



II. DILUTE BOSE GAS IN A HARMONIC TRAP 

The basic physics of a uniform dilute Bose gas has been well understood for over 30 years |p|JTo[. The recent 
intense renewed interest in such systems arises from the remarkable experimental demonstration of Bosc-Einstein 
condensation of dilute ultra-cold alkali atoms in harmonic traps p^l6|. The presence of a trap alters the physics 
in several important ways, and it is essential to reformulate the theory to incorporate the many qualitatively new 
features Q . 



A. Ideal Bose gas in a harmonic trap 

In the usual case, the trap is an axisymmetric harmonic-oscillator potential 

V tI (r) = V tl (r X , z) = \M(u\r\ + lo 2 z z 2 ), 



(2.1) 



where the coordinate vector r is expressed in cylindrical polar coordinates (r±, <j>, z). Here, u>± and oj z are the radial 
and axial angular frequencies, and it is common to introduce the "anisotropy parameter" 



A 



LU±_ 



A particle with mass M in such a potential has a Gaussian ground-state wave function 



ij)g{rx,z) oc exp 



z 
d 2 



with radial and axial dimensions 



Mu)\ 



and d, 



Mu) z 



A 1 / 2 ' 



(2.2) 



(2.3) 



(2.4) 



which are the relevant oscillator lengths (the mere presence of the trap introduces a new characteristic length scale 
that will play a significant role in the physics of these systems). The anisotropy in the ideal-gas density profile is 
simply 



dl 1 e 

—j- = — tor a non- interacting gas, 



(2.5) 



where A 3> 1 represents a "pancake" and A -c 1 represents a "cigar." 

What is the transition temperature T c for the onset of Bose-Einstein condensation in an ideal gas in such a trap? 
For many qualitative discussions, it is convenient to define a mean frequency uiq — (w^Wj) 1 / 3 and a mean oscillator 
length do = (d^dz) 1 ^ 3 . Consider the classical limit, with fc^T ^> Tiloq. Instead of the ground-state density |i/> ff (r)| 2 , 
the classical density profile has the Maxwell-Boltzmann form 



i(r) oc exp [-/3Vtr(r)] = exp 



Vtr(r) 



(2.6) 
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where (3 1 = k B T. For simplicity, consider a spherical trap (A = 1), when the density 



n(r) ocexp ( --2- ) , (2.7) 
R T/ 



is isotropic with a mean thermal radius Rt given by 

Br " M^- = d °^7 >>d °- (2 - 8) 



In this classical limit, the density is again Gaussian, but its mean dimension Rj> necessarily exceeds the mean ground- 
state radius do- Correspondingly, the mean number density is of order n ~ N/R T . 



Recall the estimate k B T c ~ h 2 n 2 ^ 3 /M in eq. (1.2) for the onset of quantum degeneracy in an ideal Bose gas. A 



combination of these expressions readily shows that 

k B T c ~ N^hujo, (2.9) 

and a detailed calculation gives the more accurate value k B T c = [((S)]" 1 / 3 ^ 1 / 3 ^ » 0.941iV 1 / 3 ?iwo- Since N typically 
is large in most cases of interest, the Bose-Einstein transition indeed occurs in the classical regime (k B T c ^> TvjJq). 
For T < T c , a Bose condensate forms in the lowest single-particle state of the harmonic oscillator potential, whose 
width is the oscillator length do, much smaller than the thermal width Rt of the non-condensate for T < T c . The 
appearance of this sharp narrow spike in the particle density served as a dramatic signal of the first Bose condensation 
in a dilute alkali gas |I3 . 

In a typical situation (see, for example, ref. |l8|), the trap has a frequency uo/2ir ~ 100 Hz and an oscillator length 
do ~ 2 /xm, containing N ~ 10 6 particles. The ground-state energy corresponds to a temperature hujQ jk B ~ 5 X 10~ 9 K. 
In this case, 7V 1/3 ~ 100, so that T c ~ 0.5 x 10~ 6 K; at this temperature, the ratio R T /d Q - TV 1 / 6 is of order 10. 



B. Effect of interactions on a trapped Bose condensate 



As in the case of a uniform condensate, the Bose gas is assumed to be dilute, in the sense that the interparticle 
spacing I ~ n^ 1 / 3 is large compared to the s-wave scattering length a and the range of the interactions. Th us th e 
interparticle potential is again approximated as a pseudopotential V(r) « g5^(r), with g w Airfi 2 a/M, as in eq. ( 1.22 ). 
Typically this low-density condition is well satisfied, with a ~ a few nm and mean number density n ~ 10 20 m -3 , 
implying I ~ a few xl0~ 7 m and a « I. In addition, the condition o C d holds whenever na 3 <C 1, since the 
existence of a macroscopic Bose condensate implies that I <C do- 

The first question is how to generalize the Bogoliubov theory to describe a non-uniform Bose gas in a trap. Although 
various approaches lead to essentially the same final results, the treatment here 19 emphasizes the field-theory aspects 
of the many-body hamiltonian 



H = 



d 3 r 



where 



T = — 



rv 

2M 



(2.10) 



(2.11) 



is the kinetic-energy operator, Vt T is the harmonic potential given in eq. (2.1), and the two-body potential has been 
replaced by a short-range pseudopotential V(r) ~ gS^(r). Here, ip and ft are field operators that obey bosonic 
commutation relations 



^(r),^!-') = <5 (3) (r-r'), ^(r),-0(r') = ^(r),^^') 
As an aside, the field operator for the uniform system has the simple form 

ft*. 



0. 



e lkr a k , 



(2.12) 



(2.13) 
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using the normalized plane- wave single-particle wave functions. More generally, given any complete set of normalized 
single-particle wave functions {Xj( r )}i the corresponding field operator is a sum over all normal modes 



3 

where a,j is a corresponding bosonic annihilation operator for the single-particle mode j. 

In the Bogoliubov approximation for a non-uniform system, one single mode is assumed to be macroscopically 
occupied. The exact field operator is therefore separated into two terms, 

$(r) « #(r) + 4>(r), (2.15) 



where <I> is the (large) condensate wave function [analogous to the Bogoliubov prescription in eq. (1.27)] and <j> is 
the (small) operator characterizing the remaining non-condensate. Formally, the condensate wave function \& = (ip) 
is an ensemble average of the full field operator in a particular ensemble that breaks the symmetry with respect to 
the phase [ pL0| | , and the remaining fluctuation part (f> can be defined by the condition (</>) = 0. In particular, the 
condensate density no(r) is given by 

n (r) - |*(r)| 2 , (2.16) 

and the total number of condensate particles is 

N = [ d 3 r |*| 2 . (2.17) 



It is most convenient to introduce the grand canonical ensemble at temperature T and chemical potential fx, studying 
the modified hamiltonian 

K = H - (j,N, (2.18) 

where N = J d 3 rip^ is the number operator. For a uniform system, the condensate number density uq is simply a 
spatial constant. For a non-uniform system, however, the determination of vfr and no is in general quite complicated. 



In analogy with the leading term of the Bogoliubov hamiltonian in eq. (1.30), where all the creation and annihilation 



operators were replaced by the condensate contributions y/No, the condensate part of the modified hamiltonian 

Ko = H -fiN = J d 3 r[y*(T + V tI - fi)^ + \g 1<*1<**1<] (2.19) 

is assumed to dominate the physics at low temperature. 

In the grand canonical ensemble at zero temperature, equilibrium corresponds to the minimum value of the ther- 
modynamic potential, which is the ensemble average of the modified hamiltonian ( K) . F or the present situation 



of a macroscopic non-uniform condensate, this quantity is approximately Kq in eq. (2. IE), and the only available 
parameter to vary is the condensate wave function itself. Thus, the value of the integral Kq must be stationary under 
the variation \&* — * ty* + 6 iff* |L9]]. The corresponding Euler-Lagrange equation 

(T + V tr - U )* + 5 |*| 2 * = (2.20) 



is known as the Gross-Pitaevskii (GP) equation |20j,|2l|]. This non- linear Schrodinger equation for the condensate 
wave function is expected to hold at low temperature; it is similar in form (but not in spirit) to the non-linear 
Ginzburg-Landau equation for the superconducting order parameter that applies near the superconducting transition 
temperature (see, for example, ref. Q, sec. 45). 

It is convenient to introduce the "Hartree" potential energy Vjj(v), which is the interaction potential of a particle 
at r with all the other condensed particles 

V H (r) » J d 3 r' V(r - r') n (r') « gn (r) = ,g|*(r)| 2 , (2.21) 

so that the GP equation assumes the equivalent form 

(T + Vu + V H - n) * = 0. (2.22) 
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The interaction energy of the condensate is the last term of eq. ( 2.19| ) 



E int = (V int ) = J d 3 r i 5 |*| 4 = i J d 3 rV H (r) |*(r)| 2 = \{V H ) , (2.23) 

where the angular bracket here denotes an integral (• • -)o = / d 3 r \f r * ■ ■ ■ \P over the ground-state condensate wave 
function; note that -Eint is half the mean Hartree energy of the condensate. 
In time-dependent form, the GP equation becomes 

ih— = {T + V tI + V H )*, (2.24) 
at 

where the condensate wave function is interpreted as an off-diagonal matrix element of the destruction operator 
\t(r, t) = (N — l\tp(r,t)\N). The time-dependence of the operator is given by the usual Heisenberg picture ^>(r, t) — 
cxp(iHt/h) ip(r) exp(—iHt/h). The time-dependent phase factor in the matrix element thus involves the change in the 
ground-state energy E a (N — 1) — E g {N) when a particle is removed from the system, which immediately reproduces 
the original GP eq. ( ^20|) because fj, « dE g /dN at T = 0. 



The ground-state solution of the GP equation satisfies several important rel ations . The first is the ground-state 
energy E g , which is the expectation value of the condensate hamiltonian in eq. ( 2.19| ) 



E g = (ffo)o = <T)o + (V tt )o + (Mnt)o = (T) + (V tt )o + ?(V H )o, (2-25) 

Second, the expectation value of the GP equation itself yields an expression for the chemical potential 

N of i - (T) + (V tr ) + (V H )o = (T)o + (Vtr)o + 2(M nt ) . (2.26) 

Finally, consider a dimensionless scale transformation on the coordinates r — > ^r. The scaling of the various terms in 
Ho follows by inspection to give 

E g -> C 2 (T)o + C 2 (Vtr)o + C 3 (V int ) ; (2.27) 

since this scale transformation cannot change the ground-state energy, it follows that (dE g /d() |^ = i = 0, which readily 
yields the virial theorem (22| 

2(V/ tr )o - 2(T) + 304rt)o = 2(T) + |(Vh) . (2.28) 
A combination of this identity with Eqs. ( p. 25 ) and ( p. 26 ) provides the simpler expressions 



£ fl = iCr>o + !<Vtr)o, (2.29a) 

N of i = -l(T) + l(V tI ) (2.29b) 
that do not involve the Hartree or interaction energy. 

C. Basic physics of the Gross-Pitaevskii equation 



The presence of the trap introduces a new characteristic energy hu>o and a new length scale do = tJK/Mojq. The 
repulsive interparticle interactions tend to expand the condensate, so that the true mean size Rq generally exceeds the 
size d of the ideal Bose gas in the Gaussian ground state. To estimate the dimensionless expansion ratio R = R /d Q , 
consider the order-of-magnitude of the various terms in the ground-state energy. Since the ground-state wave function 
is nodeless, the kinetic energy per particle has the order-of-magnitude 

<2 - 30 *» 



Similarly, 

(Ftr)o - hw Q R 2 , (2.30b) 



No 
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and 



(Vff)o h 2 an ti 2 aN hu; Na 

^~ 9n ~^r~MBZ~i?^> (2 - 30c) 

with the new dimensionless parameter Na/do that specifically reflects the presence of the trap. Although the ratio 

a/do is generally small, of order ~ 10~ 3 , N is often of order 10 6 -10 7 , so that the product Na/do can itself be large, 
of order 10 3 -10 4 . 

The total energy per particle is the sum of these three contributions 

E a + ( 1 9 Na 1 \ , 

TT~foo h^+^ + ^l^ - (2-31) 



N \R 2 d R 3 

with the dimensionless expansion ratio R as a variational parameter. Two limits are easy to analyze 



1. If Na/do < 1 (a nearly ideal Bose gas), the last term of eq. (2.31) can be omitted, and the minimum ground- 
state energy occurs at R = 1, which is just the familiar variational solution for the simple harmonic oscillator with 

(T) = (Vtr) . 

2. In the opposite limit that Na/do 3> 1) the repulsive interactions are crucial (even though the gas is assumed 
to remain dilute with na 3 < 1). In this case, the kinetic energy (oc R~ 2 ) can be neglected because the condensate 
is significantly larger than its ideal-gas dimension, and the minimum ground-state energy occurs when the last two 
terms in eq. ( |2 .31 ) are comparable, namely when 



R 5 = Rl/dl ~ Na/do > 1. (2.32) 

Before proceeding with an analysis of this very important scaling relation that holds in what has become known 
as t he Th omas-Fermi (TF) limit, it is helpful to recall the healing (or coherence) length £ = (8%na 

)- 1 / 2 defined m 



eq. (1.45) for a uniform unbounded condensate. Physically, £ arises from the balance between the kinetic energy 
h 2 /2M£ 2 and the Hartree energy gn = Anh 2 a/M. When the bulk condensate wave function 'J is perturbed locally 
(for example, an otherwise uniform Bose gas in a half space bounded by an external plane), the length £ characterizes 
the distance required for \& to heal back to the previous background value, justifying the name "healing length." 

The presence of a closed boundary (with an additional associated length scale) complicates the situation, even in the 
absence of an external confining potential like Vt T . One simple and exactly soluble example is a three-dimensional dilute 
Bose gas with uniform areal density N/A in the yz plane and confined to a slab \x\ < \L by rigid transverse walls p3| . 
Here the real condensate wave function obeys a one-dimensional GP equation — (h 2 /2M)ty" — + g 1 ^? 3 = 

and vanishes at x = ±\L. The nearly ideal limit can be characterized by the condition £ 2 ~ AL/Na 3> L 2 , when 
the condensate density uq(x) = \^(x)\ 2 is everywhere small and proportional to cos 2 (wx/L). In this case, the length 
scale L characterizes the spatial variation of n , so that L itself plays the role of the healing length in the nearly 
ideal limit. As N increases (and £ decreases), however, the repulsive interactions act to flatten and spread the central 
density maximum. In the opposite limit ( 2 < L 2 , the detailed solution shows that the condensate density no forms a 
plateau with the density falling to zero at the walls in a skin thickness of order £. This simple model shows that the 
physical healing length in a confined Bose condensate should be taken as the smaller of the two lengths L and £ 

A similar situation holds for a dilute three-dimensional Bose gas in a harmonic trap, apart from the additional 
complication that the condensate's volume ~ i?g can exceed the ideal-gas volume ~ dp when the interactions are 
sufficiently strong. The squared coherence length £ 2 ~ 1/na ~ Ro/Na provides an instr uctive way to think about 



the various limits discussed above in connection with the variational energy in eq. ( 2.31 ). In the nearly ideal limit 
(Na/do <C 1 and Ro ~ do), it is easy to see that £ 2 3> d 2 ,; the ground-state Gaussian condensate density varies on 
the length scale do, which thus acts as the healing length. As N increases, however, the coherence length £ shrinks, 
becoming comparable with the ideal-gas condensate dimension do when Na/do is of order 1. 

In the opposite (TF) limit Na/do ^ 1, the coherence length £ falls significantly below do, and £ itself now plays the 
role of the healing length (for example, near a hole produced by a localized blue-detuned laser beam). In addition, 
the characteristic condensate dimension Ro expands beyond do- To analyze this last result in detail, note that the 
definition of the coherenc e len gth implies that £, 2 R 2 ~ Ro/Na, and use of the preceding qualitative TF relation 



Rl/dl - Na/do horn eq. (2.32) yields 



eRl ~ 4 5 W- ~ & or, equivalently, -f ~ < 1. (2.33) 

Consequently, in the TF limit, the oscillator length do is approximately the geometric mean of the small healing length 
£ and the large condensate radius Ro- 
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If the gas remains dilute with na 3 <C 1 , the healing length is considerably larger than the interparticle spacing. The 
situation changes in the extreme large- N (high-density) limit with na 3 ~ 1, for the healing length £ then becomes 
comparable with the interparticle spacing I ~ n -1 / 3 . Although this dense regime is probably unattainable for the 
trapped alkali gases because of three-body recombination effects, it is worth mentioning that bulk superfluid 4 Hc 
does have na 3 ~ 1, for its healing length £ is comparable with an atomic dimension (and thus with the interparticle 
spacing) . 

The original Bogoliubov theory of a dilute uniform Bose gas with na 3 <C 1 has no independent length scale for the 
dimension of the condensate. In contrast, the new dimension do for the trapped Bose gas correspondingly yields two 
distinct physical regimes (although both are dilute with na 3 < 1): 

1. the nearly ideal regime with £ 3> do and Na/do <C 1; 

2. the interacting regime with £ -C do and Na/d ^> 1. 

This comparison of the uniform and trapped dilute Bose gases emphasizes the qualitatively new physics associated 
with the presence of the trap. 



D. The behavior of the condensate for large N 

In the limit Na/do ^ 1> the GP equation ( [2.2C ) can be simplified by omitting the kinetic energy p^pq ], yielding 
a simple algebraic equation 

(V tI + V H - = 0. (2.34) 

Either W = or the Hartree potential and condensate density are given by Vh = g\^\ 2 = A* — Vtt- In the present case 
of a harmonic potential, the approximate condensate density has a parabolic profile 

|, W P . „„(r> . „«,) (2.35) 

where denotes the unit positive step function. Here, 

n(0) = »- = -^r (2-36) 



is the central density in the trap and 



Rl = J^ and R 2 Z = ^- (2.37) 



MuP\_ " " z Mui 



characterize the radial and axial semiaxes of the ellipsoidal condensate density. In this large- N limit, the slowly 
varying trap potential determines the local density [often called the Thomas-Fermi (TF) limit in analogy with the 
corresponding treatment of the electron density in large atoms]. Note that the anisotropy of the condensate in this 
TF limit is now given by 



Rl uj 2 ± I 



-#- = -4 = -t in TF limit, (2.38) 
R, uji \ z 



differing from the anisotropy 1/A in eq. (2.5) for the near- ideal condensate. 

Typically, the fractional depletion N'/N = 1 — No/N is small, and it is sufficient to set No « N and hence no ~ n. 
In this case, the normalization of the condensate wave function requires 

J d 3 rn{r)K J d 3 r\^{v)\ 2 k N. (2.39) 

A straightforward calculation pa] with the TF density from eq. ([Oq) yields 



where Rq = Rj_R z - This result can be rewritten in dimensionless form as 
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which has many important implications for the behavior of the condensate in the large- TV (or TF) limit. 

1. The dimensionless expansion factor R grows relatively slowly, like TV 1 / 5 , which thus characterizes the expansion 
of the mean condensate radius Rq. Equivalently, for a given radius Rq, the total number increases like N oc Rq. 

2. The central density n(0) is of order N/Rq, and the preceding dependence implies that n(0) oc Rq. 

3. The chemical potential is fj, « ^TlujqR 2 , with R ^> 1, so that fj, far exceeds the ground-state energy of a single 
particle in the trap. 

4. The ground-state condensate energy E g w |(Vt r )o = |-^ w o( r2 )o = j^fouoR 2 ^ oc TV 7 / 5 follows from the virial 
theorem in eq. ( 2.29a[ ) and direct evaluation of (r 2 )o, since (T)q is negligible in the TF limit. 



5. The TF chemical potential given above in item 3 is consist ent wi th the general thermodynamic identity /i = 
(dEg/dN) « lEg/N oc TV 2 / 5 , as well as the virial theorem in eq. fl2.29b| ). 



6. It is natural to define the healing length for the trapped condensate in terms of the central density, with 
£ = [87rn(0)a] _1 / 2 . In the TF limit, this choice implies that 



£i?o = dl, or, equivalently, J- = ^- < 1, (2.42) 

«o Ro 



quantifying the previous qualitative TF conclusion from eq. (2.33) that the bare trap size do is the geometric mean 
of the healing length £ and the mean condensate radius Rq . 

7. This whole TF analysis is easily generalized to an arbitrary anisotropic (triaxial) trap with u x =/= u y =/= lo z by 
reinterpreting the mean frequency as ujq = us x u y u) z , with similar expressions for dy = d x d y d z and Rq = R x R y R z . 



Recall that a uniform Bose gas has the speed of sound s = J gn jM from eq. ( 1.44 ). In a trap, the central density 



n(0) provides the corresponding estimate s « y/gn(Q)/M, and item 2 above shows that the central density n(0) oc Rq 
scales with the square of the condensate radius Rq in the TF limit. Thus, the TF speed of sound in a trapped 
condensate increases linearly with Rq. With item 6 above, it is easy to find the quantitative expression 

h Rn h , 

(2.43) 



V2M dl V2M£' 

This linear dependence on the condensate radius has the following remarkable implication for the lowest-lying com- 
pressional modes of the condensate. The fundamental mode has a wavelength of order Rq, and the lowest frequency 
will therefore be of order u> ~ s/Rq ~ h/Md^ ~ ujq, namely the mean trap frequency ojq. Consequently, the low- 
est normal-mode frequencies of the condensate should be independent of the condensate radius Rq (and hence the 
condensate number N) in the TF limit of a large condensate. Several experimental studies [^6|,^7| have verified this 
prediction (sec. |l| below contains an account of the theoretical analysis 2j|). 

It is also instructive to consider the Landau critical velocity v c for the onset of dissipation in the condensate, which 
is simply the speed of sound in the Bogoliubov description of a uniform dilute Bose gas. In the present case of a 
trapped condensate, it is not obvious how to send an impurity rapidly through the sample, but a rigid rotation with 
angular velocity can accomplish much the same effect. The equatorial speed is ~ QRo, and the onset of dissipation is 
expected when this speed equals v c w s. Thus the critical rotation speed should be of order Q c ~ v c /Rq rs s/ Rq ~ ujq 
as shown above. This result means that the condensate should become unstable if the rotation speed exceeds any of 
the trap frequencies (a similar instability occurs for a particle in a harmonic potential). 

E. Effect of an attractive interaction 

Bose-Einstein condensation has been observed |l6| ] in 7 Li, which is known to have a negative scattering length a, 
corresponding to an attractive interaction. For a uniform dilute Bose gas, the corresponding energy spectrum in eq. 



(1.42) becomes imaginary at long wavelengths (see fig. 3) 



El = {jm) fc2 (fc2 ~ 16 ™ H) ' (2 - 44) 

implying an instability for those modes with k 2 < 167rn|a|. Equivalently, the speed of sound becomes imaginary. For 
a trap, however, the wavenumber cannot be arbitrarily small, and the minimum value is of order fc m j n ~ tt/Rq. Hence 
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the system can remain stable if tt 2 /Rq > 167rn|a| 
critical number 



Since the density is of order n ~ N/Rq, this analysis predicts the 



c ~ 16 lol ~ 16 lap 



(2.45) 



because i?o ~ do for such small values of the ratio A c |a|/do ~ 7r/16. For A < A c , the trapped condensate with 
attractive interactions should remain stable because of the positive kinetic energy arising from the confinement; for 
A ^ N c , however, the attractive interaction energy predominates, and the condensate should collapse. 

It is not difficult to make a variational estimate of the critical condensate number N c for a spherical trap with 
attractive interactions. As a simple model, take a Gaussian trial function [a similar variational estimate was used by 
Baym and Pethi ck |B5| in their initial study of the GP equation for repulsive interactions; compare also the qualitative 
estimate in eq. ( |2.31| )f, 



ty(r) oc exp 



with P as the variational parameter; the condensate radius is 
ground-state energy is readily evaluated to be piT 



(2.46) 



f3dp with j3 < 1 for attractive interactions. The 



1 



1 



2 \p 



2 A|a| 1 
7T do~W 



(2.47) 



and fig. 4 shows this function of /? for several values of the interaction parameter N\a\/da. Evidently, E g (/3) has no 
global minimum, for it becomes arbitrarily negative as (3 — > 0. Nevertheless, the energy does have a local minimum if 
A < N c . A straightforward analysis shows that the minimum disappears at the critical value N c \a\/do ~ 0.671, and 
that the critical condensate radius is reduced by a factor (3 C = w 0.669. For comparison, a numerical study of 

the GP equation |32| yields the value N c \a\/d n w 0.575, which differs from the variational estimate by sa 17%. 



F. Surface region of Thomas-Fermi condensate 



For a spherical trap, the slope of the TF condensate wave function 



* T p(r)cx Jl-^Q(R -r) (2.48) 

becomes infinite at the surface (r — > Rq), leading to logarithmic singularity in the kinetic energy. Consequently, 
it is essential to introduce a boundary layer near the surface, where the true wave function is modified from its 
approximate TF form |p3| . To study this surface region, it is convenient to focus on the region \r — Ro\ <C -Ro, writing 
r = Ro(l + XS), where A is a new scaled dimensionless variable measured from the TF surface, and 5 <C 1 is the 
dimensionless thickness of the boundary layer, measured in units of Rq. Equivalently, A = (r — Ro)/ (Rod). With this 
change of variable, the GP equation for a spherical trap 

h 2 d 2 T h 2 d T / 1 , , 9 9 \ 4irh 2 a , T , T /„ ,„s 

* + ^Mw r 2 - u) *H * 2 * = (2.49) 

2M dr 2 Mr dr ^ 2 ' M 1 ' 



takes the approximate form 



!r ' "' 2 * + Mo;^X* + ^^|*| 2 ^ = 0. (2.50) 



2MR 2 S 2 dX 2 ''-oo—' M 
The balance between the first two terms gives the relation 



^m^pp!^ 4 ^ 4 ' 5 ' (2 - 51) 



verifying that the boundary layer is indeed thin and showing that 5 oc R 4 / 3 scales with the condensate number 
like A -4 / 15 . In contrast, the balance between the last two terms shows that the condensate density l^j 2 near the 
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surface is small (of order S). A detailed study [|33|-p5| of the condensate wave function in the surface region involves 
a non- linear differential equation. Its solution matches the TF behavior in the interior of the condensate for X <C — 1 
and vanishes exponentially for X ^> 1, thus providing a uniform extension across the boundary layer. The resulting 
kinetic energy per particle 



(T) 
N 



5 h z 



2 MRl 



hi 



Ro 
1.3d Q 



(2.52) 



exhibits the role of the boundary layer in cutting off the logarithmic singularity. This contribution to the total energy 
is much smaller than the other dominant contributions (Vt T )o/N j^MljqRq and (Vi n t)/N = |{Vj?)/iV « IMlJqRq, 
whose sum gives the previous expression E g /N sa j^Mlo^Rq. 



G. Excited states of a trapped Bose condensate 



In the case of the uniform stationary dilute Bose gas, the condensate wave function is a spatial constant. This simple 
structure allows a complete explicit solution for the elementary excitations, which are the Bogoliubov quasiparticles. 
For a non-uniform trapped Bose condensate, however, the situation is more complicated because the determination 
of the condensate wave function \&(r) already constitutes a significant problem. Nevertheless, given 'J', a well-defined 
procedure leads to the corresponding equations for the small-amplitude normal modes of the condensate. These 
equations were first studied by Pitaevskii [pl| in the context of a long straight vortex line in an unbounded condensate, 
and they have sub sequ ently been generalized as operator equations in ref. |l9| . 

Recall that eq. ( [2.15] ) separated the field operator into a (large) condensate part and a (small) fluctuation operator 



v> = * + 



(2.53) 



where "J = (ip) and (</>) = 0, and the angular brackets denote an ensemble average in a restricted grand-canonical 
ensemble that breaks the symmetry with respect to the phase |l(|. The derivation of the GP equation for the 
condensate wave function retained only the domi nant terms involving 'J' , giving the leading condensate approximation 
Ko to the thermodynamic potential in eq. ( 2. IS ), thereby ignoring the contribution of the fluctuation operators. To 
include the first corrections, it is necessary to retain the fluctuation contributions <j> and (jy to the original field 
operators; they obey the approximate bosonic commutation relations 19 



r),«£V)l «<5C3)( r -r'), U(r),0(r')l = U+(r), ^( r ')l = 



(2.54) 



Substitution into the full modified hamiltonian K = H — fiN from eqs. (2.10) and ( 2.18| ) gives the previous Kq plus 
a quadratic correction 



K' 



d 3 rft (T + V tI + 



2irh 2 a 
M 



d A r 



4|*| 2 t + # 2 <^t + (^*) 2 



(2.55) 



The linear contribution to K vanishes identically because K is stationary with respect to small changes around the 
GP solution 4"; this condition merely requires that 4' obey the appropriate Euler-Lagrange equation. 

The mean number density in the ground state is given by n(r) = (^M '(r)ip(r)) . The expansion in eq. (2.53) shows 
that n(r) = ^(r)! 2 + (<fi^ (r)<fi(r)) , because, by assumption, the linear contributions have zero ground-state ensemble 
average. Consequently, the mean total number of particles is a sum of condensate and noncondensate contributions 
N = N + N', where N = J d 3 rn (r) and 



N' = / d 3 r(ft(r)j)(r)). 



The dynamical equations for the fluctuation operators follow immediately from the commutators with K' : 



and ih ■ 



dt 



(2.56) 



(2.57) 



A straightforward calculation yields the coupled operator equations 
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ih = ( T + Vtr - n- 



8irTi 2 a 
M 



1*1 



M 



2 It 



(2.58a) 



1*1 



47r^ 2 a 
M 



(**) 2 0. 



(2.58b) 



In analogy with Bogoliubov's treatment of a uniform dilute Bose gas, assume that the fluctuation field operator has 
an expansion in a set of self-consistent normal modes 



M)=]T' [ U] {v)a 3 {t) - v*{v)a]{t) 



(2.59) 



where the primed sum runs over all the excited states. Here, the quasiparticle operators have harmonic time depen- 
dence (Xj(t) — ctj exp(—iEj /%) and ctj{t) = o^j exp(iEj /U), obeying bosonic commutation relations at equal times. 
Substitution into eqs. (2.58) leads to the coupled "Bogoliubov equations" pq] 



T + ~ H ■ 



iirh 2 a , 



M 



1*1 



4:irh 2 a 
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(2.60a) 



T+Vt t -ii- 
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(2.60b) 



that determine the pair of quantum eigenamplitudes (in effect, they are wave functions) Uj(r) and Vj(r) along with 
the associated energy eigenvalue Ej for the jth eigenstate. 

This set of equations can be rewritten in matrix form as MLij = Eja z Uj, where A4 is a 2 x 2 matrix differential 
operator, Uj is a two-component vector with elements Uj and Vj, and a z (the usual Pauli matrix) acts as a metric 
tensor in the scalar product / d 3 rU\a z Uj. The diagonal elements M\\ — M22 = T + V tI — fi + 2Vh of the matrix 
M consist of hermitian operators (in the usual quantum-mechanical sense), and the off-diagonal elements M21 — 
M.\2 = — {Airfi 2 a / M)^> 2 are complex-conjugate functions, so that = M.. With suitable boundary conditions, 
an integration by parts shows that J d^rU^MUj = j d 3 r {MlAkfi Uj- As in the usual Sturm-Liouville theory, it 
follows |l9| that the eigenvalues Ej are real and that eigenfunctions for distinct positive energies Ek ^ Ej are 
"orthogonal" with J d 3 r {u* k Uj — v^Vj) — 0. For Ej > [see eq. (2.62) below], the proper choice of normalization 



d 3 i 



(h 



= 1 



(2.61) 



1 in eq. (1.32) for the uniform system; it ensures that the quasiparticle 



is analogous to the condition u 2 — v\ 
operators obey the usual bosonic commutation relations. Although the Bogoliubov equations describe the excited 
states, the condensate wave function ^ itself acts as a "zero- mode" solution |L9|]; its conjugate operator completes 
the set of states and is important in the quantum dynamics of the phase of the condensate P3,pq] • 



In addition, for every solution Uj 
-Ej and normalization J d 3 r [\vj 12 



and Vj with energy Ej, there exists another solution v* and u* with energy 



; 3 | 2 ) = — 1; this second set corresponds to interchanging at and otj, in effect 
exchanging quasiparticles and quasiholes. Although the Bogoliubov equations describe bosonic amplitudes, their 
structure is very reminiscent of the Dirac equation, where the coupled amplitudes again include particles and holes 
that have positive and negative energies |39| (the qualitative differences appear only when each is reinterpreted as a 
quantum field theory with distinct bosonic or fermionic commutation relations). 

Equation ( [2.56 ) shows that the total non-condensate number AT' is the ground-state expectation value f dP rffify. 
Since the quasiparticle operators represent non-interacting bosons, use of the normal- mode expansion (2.59) yields 
N' = Ylj'Nj, where iVj = / d 3 r\vj\ 2 is the ground-state occupation of the jth eigenstate. In addition, direct 

substitution of eq. (2.5E) into eq. ( |2.55 ) for the quasiparticle modified hamiltonian K' and use of eq. (2.60) gives the 
simple and intuitive expression 



k' = 



EjN> } 



E E 3 a j a 3i 



(2.62) 
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which is a direct generalization of the Bogoliubov quasiparticle hamiltonian (1.39) for a uniform system. The first 
term is the ground-state thermodynamic potential K g = (K') = — J2j'Ej iVj, and the second shows that the excited 

states consist of various numbers of non-interacting bosons with energies Ej and occupation- number operators ctjCtj- 
Since these operators have non-negative integers as their eigenv alues, stability of the ground state requires that Ej > 



for all eigenstates j that obey the proper normalization in eq. (2.61). 

At low but finite temperature in the grand canonical ensemble, the appropriate (unnormalized) weight factor for the 
excited states is exp(— /3K'); the average of any operator O is simply the trace (O) = Tr[0 exp(— (3K')]/Ti[exp(— (3K 1 )] 
with the same weight factor (normalized by the grand partition function Z — Tr[exp(— (3K')\). For example, the non- 
condensate density is 

»'M = E' {/(^)Mr)| 2 + [1 + m)] Mr)| 2 } , (2.63) 

3 

where f(Ej) = [exp((3Ej) — 1] 1 is the Bose- Einstein distribution function for the jth excited state. At zero temper- 
ature, the spatial integral gives the preceding result for the total non-condensate number. 

It is instructive to specialize the present Bogoliubov quasiparticles to a uniform system. In sec. |[ the condensate 
number Ao was eliminated in favor of the total number N, leading to a number-conserving description (with no 
chemical potential \x in the hamiltonian); in contrast, the theory here relies on the grand canonical ensemble at fixed 
temperature T and chemical potential fx, leading to a modified hamiltonian K = H — /iA. In this latter formalism, 
the ground-state thermodynamic potential per unit volume for the uniform condensate is Kq/V = \gn\ — fino, 
where no = No/V is the condensate density (the kinetic energy vanishes for a uniform condensate). In equilibrium 
at fixed chemical potential and temperature (here, T = 0), the system adjusts all free parameters to minimize the 
thermodynamic potential. In the present case where the non-condensate is neglected, only the condensate density no 
remains undetermined, and the condition (dKo/dno)^ = immediately gives the expected result u = gno, which is 
the Hartree energy for adding one particle. The corresponding solutions of the Bogoliubov equations are plane waves 
Uk(r) — u k exp(ik • r) and i>k(r) = v k exp(zk • r). Use of the preceding expression for u in the Bogoliubov equations 
( [2.60 ) readily reproduces all of previous Bogoliubov results, for example, E k = ^2gn e k + (e k ) 2 



More generally, the Bogoliubov equations provide a flexible formalism to study the effect of a non- uniform condensate 
wave function. A particularly interesting example is a uniform system in which the condensate moves with velocity 
v = hq/Al (discussed briefly at the end of sec. |). The condensate wave function has the form \& cx e !q r , and the 
GP equation then shows that the chemical potential becomes /i = e° + gno = ^Mv 2 + g n o- The structure of the two 
Bogoliubov eqs. ( [2.60 ) implies that the pair of amplitudes has the form 



u(r)\ _ f u k e^ r e jk r 
v(r) ) ~ \v k e- tcl r e lk -' 



(2.64) 



and the resulting eigenvalue is given by E^ — ?ik ■ v ± E®, where E® is the Bogoliubov energy for the stationary 
condensate. Substitution back into the Bogoliubov equations shows that the ± sign here corresponds to the normal- 
ization uf, — v\ = ±1, so that the bosonic commutation relations require the choice E^ = fi,k • v + E^ (as also follows 
by requiring continuity with the small- v limit). For a detailed discussion of these solutions, along with the associated 
mass current and superfluid density (as well as the generalization to a mixture of two distinct interacting species), 
see, for example, refs. |[l^ , ^0| . 

Consider an anisotropic harmonic trap with arbitrary frequencies uj a (a — x,y,z), and let ^ denote any exact 
solution of the GP equation. In general, the Bogoliubov equations are difficult to solve, but they do have three exact 
solutions (the "dipole" modes) in which this general condensate oscillates rigidly in each orthogonal direction with 
frequency uj a . To construct these explicit solutions [ffl| , define the familiar raising and lowering operators 



t _ 1 (x a Pa\ _ J_ / Xq_ _ , d \ 



(2.65a) 



It is not difficult to prove that the corresponding dipole solutions have the explicit form 

(2.66) 



u a (r) \ = ( ot*(r) 
v a (r) J U„r(r) 
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For the non-interacting Bose gas, the solution of the GP equation is just the Gaussian ground state \& s = vNi^ g of 
the trap given in eq. ( [2.3] ); the associated dipole excitations are then characterized by u a oc a^ipg (namely, the first 
excited states of the oscillator) and v a = (because a n a nnihilates the ground state tp g ). 

It is worth noting that the Bogoliubov equations ( 2.6C ) have an equivalent derivation |fl2|-f44| based directly on the 
linear response of the time-dependent GP equation (although this approach leads more directly to the final equations, 
it relies only on wave functions and does not treat the field operators explicitly). Add a weak sinusoidal perturbation 
to the GP equation 



in ^ = (T + V t , + g\*\ 2 ) * + (f+e- iut + /_e Mt ) *, 



(2.67) 



and seek a solution in the form 



[#(r) + u(r)e 



■«(r)e* 



(2.68) 



Linearizing in the small amplitudes u and v yields an inhomogeneous form of the Bogoliubov equations with driving 
terms proportional to f±{r) 1 $>. A normal- mode expansion in terms of eigenfunctions uj and Vj gives precisely the 
same coupled equations as in eq. (2.6C), identifying the frequencies ujj — Ej/% as resonances in the linear response. 



III. HYDRODYNAMIC DESCRIPTION OF A DILUTE TRAPPED BOSE GAS 



Although the Bogoliubov equations, in principle, characterize all the excited states of the condensate, they are 
difficult to solve and their physical meaning is not always transparent. Thus, it is important to consider alterna- 
tive formulations, and the present section concentrates on the hydrodynamic approach that emphasizes the familiar 
concepts of density and current fluctuations. 



A. Uniform dilute Bose gas 



It is instructive to review the theory of density fluctuations in a uniform Bose gas, where the field operator has the 
familiar form 

v3(r) = £^=e* k - r a k . 

The corresponding density operator n(r) = ^ (r)ijj(v) has a Fourier expansion 

n(r)=^e^ r n q , (3.2) 
q 

where 



E 

k 



a k-q fl k (3-3) 



is an operator that removes one quantum of density fluctuation with wave number q. It is easy to see that the 
associated creation operator is given by n q = n_ q . For an ideal Fermi gas at T = 0, the operator fi q in eq. ( |3.3| ) 
removes a particle with wave number k (which must lie inside the Fermi sea with |k| < fop) and simultaneously 
creates a particle with wave number k— q (which must lie outside the Fermi sea with |k — q| > hp); the final operator 
fi q is a sum over all possible k consistent with the above restrictions. Thus the density-fluctuation operator n q is 
really a particle-hole operator that conserves particle number, and a similar characterization applies to both a general 
interacting Fermi gas and a non-condensed Bose gas. Typically, density fluctuations in a Fermi gas have a well-defined 
dispersion relation Uk, and they are often called "collective modes" to distinguish them from the "quasiparticles" that 
are associated with single-particle excitations (those resulting from the removal or addition of one particle). 

The situation is qualitatively different for a condensed Bose gas, because the condensate itself plays a special role. 
For simplicity, consider a dilute uniform Bose gas with stationary condensate at low temperature. In this case, the 



two terms arising from the Bogoliubov prescription dominate the general sum in eq. (3.3), leading to the simple 
density-fluctuation creation operator 
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m v^Vo («, ~ w,) (a* + a-q) (3-4) 

that is a coherent linear combination of a quasiparticle and a quasihole. When this operator acts on the Bogoliubov 
ground state |$), it gives 

^ q l < f> = V / ^oK-^)« q l*> (3-5) 

because, by definition, a_ q |<I>) vanishes. Apart from an overall factor, the action of the density-fluctuation operator 
n q on the Bogoliubov ground state is the same as that of the quasiparticle operator aj. . Consequently, the excitation 
spectrum associated with a density fluctuation in a dilute Bose gas is necessarily identical with that for a quasiparticle; 
as discussed above, this situation is completely different from that for a dilute Fermi gas. 
The current-density operator has a similar expansion 

j(r) = ^ [^(r)W(r) - (V^(r)) = £ j,^*, (3-6) 

with 



q 



j q = ^E( k -^)4- q «k. (3.7) 



M 

k 



In the Bogoliubov approximation, two terms again predominate with 

k q 



~ |^V^o(w 9 + «?) ("q-a-q) ; (3-8) 



as expected for a dilute gas, this operator is manifestly longitudinal (namely, along q). 
For time-dependent (Heisenberg) operators, conservation of particle number requires that 



0n q (t) 
dt 



+ *q-j q (*) = o. (3.9) 



It is clear from eq. (3.4) that the density- fluctuation operator n q (t) in a dilute Bose gas oscillates at the Bogoliubov 
frequency uj q = E q /%. 

In a bulk uniform system, the ensemble average of n q n q is called the "static structure function," with the precise 
defin ition S(q) — N^ 1 (n q n q ) . For a dilute Bose gas at low temperature, this quantity is readily evaluated with eq. 
(3.4) to give 

e a . , ( E a 



S(q) ~ (u q — v q ) (a^atL + aj_a q ) = — coth ( q ] , in the Bogoliubov approximation, (3.10) 

E q \2kbT J 

since off-diagonal ensemble averages like (a q a_ q ) vanish. At T = 0, the thermal factor is simply 1, giving the 
zero-temperature limits 

S(„_ 4 »/ "*/*•". J-**'- (3.1!) 
Eq I 1, for q£ > 1, 

where s is the speed of sound. More generally, the Bogoliubov approximation at T — reproduces the form of 
Feynman's variational approximation for the zero-temperature density excitation spectrum of supcrfluid 4 Hc Eq] 



k 2 q 2 
2MS(q) ' 



E 1 « (3-12) 



with S(q) the measured structure function. For low but finite temperature, in contrast, the static structure function 
for a dilute Bose gas has the limiting long-wavelength form 
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B. Sum rules for a uniform dilute Bose gas 



For both Fermi and Bose systems, sum rules provide a powerful and exact approach to the description of collective 
modes p6| , ^7| . Consider the response of a uniform system to a weak scalar field <fr(r,i) that couples to the density 
with a perturbation hamiltonian (here expressed in Fourier components) 



(3.14) 



The general theory of linear response shows that the induced change in the density occurs at the same wave vector q 
and frequency u>, with the explicit form 



S(n(q,u)} = -x(q,w)$(q,a;). 



(3.15) 



Here, the coefficient is known as the "dynamic susceptibility," and the explicit minus sign is added here for later 
convenience; x is really the Fourier transform of a density-density correlation function. Linear-response theory ]46| ] 
yields the explicit zero-temperature form 



x(q,w) 



j¥0 
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uj — ujjo + irj uj + (jjo + it) 



(3.16) 



where the limit 77 — ^ + is taken at the end of the analysis. Here, \j ) is an exact eigenstate of the interacting but 
unperturbed many-particle system, and ujjQ = (Ej — E^)/Ti is the exact excitation frequency of this state relative to 
the ground state. 

It is clear by inspection that x(q, ui) has poles in the complex u> plane at the points u> = ±u>jo — irj. Use of the 
formal identity 



lim 



1 



1 



iir8(x — a), 



r/^o+ x — a + irj x — 

where V denotes the Cauchy principal value of an integral, yields the explicit expression for the imaginary part 



(3.17) 



Im X (q,w) = ^^|(i|f^|0)| 2 [^-^- )-^ + ^ )]. 



(3.18) 



Note that I mx(q , to) is an odd function of u> and therefore vanishes at w = 0. Furthermore, a combination of eqs. 
(3.16) and ( 3.18| ) yields the spectral representation for the susceptibility (the response function) at wave vector q and 
frequency lu 



X(q,w) = / ; 

In particular, the limit u> — > is simply the static susceptibility 

[°° diJ> Imx(q,^) 

Xq = x(q, 0) = / ; — 



(3.19) 



(3.20) 



that characterizes the response to a static plane- wave distortion with wave vector q [the convergence factor irj is now 
unnecessary because Imx(q, 0) =0]. 

This result can be generalized by defining the ith moment of Imx 



/OO T I 
— j im X ( q , ui) = - 1 u K 10 ) i 2 K- y [1 - i-iy] ■, 



^From this perspective, the preceding eq. (3.20) is the i = — 1 moment 

M-i = - — lmx(q,w)=x q - 
71" Jo ^ 



(3.21) 



(3.22) 
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In the limit q — > 0, this quantity is related to the bulk thermodynamic compressibility because x q is the response to 
a static plane wave cx e' q r ; a detailed analysis J46| shows that 



giving the limit of M-i as q — > (this result is known as the "compressibility sum rule"). 
It is clear that the even moments vanish. The next odd moment 



(3.23) 



2 f°° 

Mi = — dw ui Imx(q, u) 

I" JO 
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j¥0 



can be related |4(| to the ground-state average of a double commutator (0 

follows from the dynamics 
leading to the final result 



(3.24) 



10). The inner commutator 



follows from the dynamical eq. (3.9) that expresses the conservation of particles since [H, n_] = iftdnl/dt — — q ■ j q , 



Mi 



(3.25) 



which is known as the "/-sum rule" 46 1. 

It is also valuable to consider the third moment 



M 3 = - 



2 f 

— / dio u 3 Imx(q, uj) 
k Jo 
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(3.26) 



which can be related to the more complicated matrix element 

(0| [[h^H], [H, [H,n^]\] |0) cx (0 | [q • j q , [H, q • j q t]] |0). 

For a dilute uniform Bose gas, a detailed analysis [ p0[ that neglects a small kinetic-energy part of order V na 3 yields 

^3 = ^ F e°(e° + 2n ff ) = ^-K 2 ! (3.27) 



where E q is seen to be the Bogoliubov energy spectrum in eq. ( |1.42| ). 

These various sum rules can serve to characterize the detailed structure of lm\(q, ui) |48|-po|| for a dilute uniform 
Bose gas. Assume that an excitation with wave vector q constitutes a single long-lived collective mode with frequency 
Uq. In this case, the spectral weight has only one 5-function, with Inr^q, u) — AqS(uj — uj q ) for positive uj. Direct 
substitution shows that 



Tl 1 2 2 A 

Mi = — = -A^ ) 



(3.28a) 



-J—F 2 

%M q 



(3.28b) 



The ratio of these results shows that w q is precisely the Bogoliubov frequency E q /%, with the weig ht fac tor A q = 
iT€ q / E q . In addition, substitution of the single-mode approximation into the compressibility sum rule ( 3.2C ) gives the 
static susceptibility 



kg 2 

Mu>\ 



Ms 2 



as q — + 0, 



(3.29) 



in agreement with the general result in eq. ( |3.23 ). 
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Finite temperature requires an ensemble average, weighted with the Gibbs factor e @ E * for the general initial state 
\i), along with the partition function Z = £\ e _/3Ej . The "dynamical structure factor" 5(q, ui) is defined as 

S(q, w)=^|(/ |^|i>| a 6 (w - u> fi ) , (3.30) 

fi 

where uifi — [Ef — Ei)/% can now be both positive and negative. For liquid 4 He, neutron scattering can measure 
5(q, (J) directly. In the Bogoliubov approximation for a dilute uniform Bose gas, it is not difficult to show that 

S(q, u)*§-{[l + ftE 9 )]8(w-^+f{E q )8(w + ^y (3.31) 

where f(E q ) is the Bosc-Einstein distribution function. The first term (representing emission of a quasiparticle with 
energy E q ) is known as a "Stokes" process; it persists even at zero temperature. The second term (representing the 
absorption of a quasiparticle with energy E q ) is an "anti-Stokes process" that requires the presence of a previously 
excited quasiparticle and thus vanishes as T — > 0. On general grounds [ ff6|| , the dynamical structure factor satisfies 
several important relations; for example, its zeroth moment is the static structure function S(q), which here gives 

S(q) = / dw S(q, w) = -f- coth — «- , (3.32) 



E„ \2k B T 



in agreement with the previous eq. (3.1C); furthermore, the first moment is just the /-sum rule at finite temperature 



hq 2 

dujujS(q,uj) = — ; (3.33) 

finally, S(q,u>) obeys a detailed-balance condition S(q, —ui) = e~P huJ S(q,ui). These relations provide non-trivial checks 
on the Bogoliubov approximation to S(q,u>) in eq. ( [3.31 ). 



C. Hydrodynamic description of non-uniform dilute Bose gas 

The condensate wave function ^ can be rewritten as 

*=|*|e iS , (3.34) 



where the absolute value is related to the condensate density |^(r)| = y/no(r). Furthermore, the condensate particle 
current 



identifies the condensate velocity as 



v (r) = ^VS(r). (3.36) 



The local superfluid velocity vo is proportional to the gradient of the phase of the condensate wave function, implying 
that V x v = 0. Hence the condensate flow is longitudinal (or irrotational) , as proposed by Landau (T^| for the 
superfluid velocity v s in H e. 



Substitution of eq. ( 3.34 ) into the GP equation ( 2.2C ) provides an intuitive hydrodynamic picture of the condensate. 



The real part gives a generalized Bernoulli's equation |20| , pi[ 

! . r 9 AirTi 2 a h 2 , ,„ 

^ = V tr + \Mv 2 + — — n a -— — V 2 ^, 3.37 

where the last term is a quantum pressure associated with spatial variations in the condensate amplitude (and 
thus the condensate density); the imaginary part is the static continuity equation for the condensate 

V- (n v ) = V-jo = 0. (3.38) 
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What determines the dynamics of small time-dependent perturbations about the static solution \l/(r)? Two equiv- 
alent approaches yield the same final results. 

1. Write the time-depende nt con densate wave function as ^(r, t) = y/n(r, t) exp [iS(r, t)] and substitute int o the 
time- depe ndent GP equation ( 2.24 ). The real and imaginary parts yield time-dependent generalizations of eqs. ( 3.37 ) 
and (3.38), which can then be linearized around the static condensate contributions |2£j] . 

2. Retain the operators <j) and cp 1 in the expansion of the field operator rp — ^ + $ from e q. (|2.53| ). This latter 
method allows a transparent treatment of the operator aspects of the hydrodynamic variables |pl| and will be used 
here. 

In a dilute Bose gas at low temperature, the total number-density operator 



naturally separates into a large condensate part no = \^\ 2 and a small non-condensate part 

ti PS + 



(3.39) 



(3.40) 



that characterizes the local time-dependent particle-density fluctuations. Similarly, the local time-dependent fluctua- 
tions in the current-density operator have the form 



n v 



n V3>', 



where 



2Mi 



(V 

*l 2 V 



(3.41) 



(3.42) 



is the perturbation in the velocity potential operator (this quantity is effectively the perturbation in the phase 
operator). Since bot h n' and $' are linear functions of 4> and <t>\ it follows that (ft 1 ) and (<&') vanish in the ensemble 
discussed below eq. (2.53), although correlation functions like (n'(r, t)n'(r' , t')} are finite. Note that n' is distinct from 
the operator <jfl<j> whose ensemble average is the non-condensate density in eq. ( 2.56| ). 

It is straightforward to verify from the equations of motion (2.58) for cp and <fy that these hydrodynamic operators 
n' and obey the linearized conservation law 



TST +V.j< = f + V.(»'v„ + „„V*')=0 



and the linearized Bernoulli's equation 



~dt 



+ v • V$' 



n'- 



4M 2 n Q 



TIqW 







(3.43) 



(3.44) 



that is a direct quantum generalization of Bernoulli's equation for irrotational compressible flow [pl| at constant 
entropy. Here, the first two terms constitute the total (or hydrodynamic) derivat ive ( see, for example, ref. ||, sec. 
48), and the last term is again the quantum (kinetic-energy) contribution. This eq. ( |3.44| ) has one significant advantage 
over the original field equations (2.58), for the chemical potential no longer appears explicitly (jj, does, of course, affect 
the condensate wave function and condensate density). 

These coupled operator equations can be rewritten with the normal-mode expansions from eq. ( p. 59 ). In this way, 
the hydrodynamic operators themselves have normal-mode expansions 



i* t 
" •i a j 



(3.45a) 



*' = £' (*">, • 'K<>, 



(3.45b) 



where the amplitudes n'j and $^ obey the same coupled eqs. ( 3.43] ) and (3.44) as the operators. Comparison with the 
Bogoliubov eqs. (2.60) shows that 



n'j = '&*Uj 



.i i 



(3.46a) 
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^ = 2^^ VS + ^' (3 - 46b) 

showing that the Bogoliubov amplitudes Uj and Vj are linearly related to the hydrodynamic amplitudes n'j and $^ 
pl|,p2[ . Thus, any solution for n'j and 4>^- gives a corresponding solution Uj and ty, and vice versa. 
For example, eq. ( 2.66| ) shows that the Bogoliubov equations have exact dipole-mode solutions u a = a\^> and 



Va = o, a ^f*, in which the condensate oscillates rigidly with frequency u> a along the three orthogonal directions (this 
holds for any solution \& of the GP equation) . Use of the explicit form of the raising and lowering operators a' a and 
a a in eqs. ( 2.65| ) and the linear relations in eqs. (3.46) readily yields the corresponding density perturbation 



n <x = 7S^— . ( 3 - 47 ) 



y/2 dx a 

which is a rigid shift of the condensate density |\I/| 2 = no, and 

,= / h f x a , dS \ ,„ 

K = -?= — -r - Mas— ■ 3.48 

The last term involving the phase S is absent if \& is real, and the velocity perturbation is then a spatial constant, 
moving 90° out of phase with the density oscillation. More generally, the last term reflects the change in the local 
velocity arising from the rigid motion of the condensate velocity. 

For a nearly ideal Bose gas, the Bogoliubov equations for Uj and Vj often provide a simpler description than 
the hydrodynamic variables; in contrast, the hydrodynamic amplitudes n'j and $^ are usually preferable in the 
large-condensate (TF) limit, where Stringari (see below) has constructed explicit solutions for the low- lying normal 
modes [P8| . As an application of the linear relation between the two descriptions, these TF hydrodynamic solutions 
immediately give the corresponding Bogoliubov amplitudes Uj and Vj in the TF limit. The spatial integral of \vj\ 2 
is the zero-temperature occupation number iVj(O) of the jth mode, and a detailed evaluation for a large spherical 
condensate shows that they are large (of order i^/dp) jyj. 

D. Hydrodynamics in the Thomas-Fermi limit 

Assume a large stationary condensate with Vo = and condensate density no = (M/4Trh 2 a)(fj, — Vt T )@(ii — Vt r ). In 
this TF limit, it is consistent to omit the quantum contribution in Bernoulli's equation, in which case the hydrodynamic 
amplitudes n'Ar)e~' lu ' it and <Fj(r)e -lUjt obey two coupled linear equations 

- iujn'j + V ■ (n V$j) = 0, (3.49a) 
, 4irh 2 a , 

- + -j^-n'j = 0. (3.49b) 

They can be combined into a single second-order equation for the density perturbation p£| 

1 



2„/ 



r " 3 ' M 

that has numerous explicit solutions. 



; + ttV • [(/x - V u ) Vn;] = (3.50) 



1. spherical trap 

In the TF limit, the chemical potential becomes \x « ^TiwqRq/ = ^AIluqR 2 , so that eq. ( |3.5C ) becomes 

uj 2 n'j + |wgV • [(Rl - r 2 ) Vn'j] = 0. (3.51) 

The spherical symmetry allows solutions of the form n' nl (r)Yi m (0,<p), where Yi m is a spherical harmonic and n[ l} is 
the corresponding nth radial eigenfunction. Solutions with different m are degenerate, and the TF frequencies |2q] 
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u'li = [I + n (2n + 21 + 3)] , 



(3.52) 



can be compared with the frequencies lu ti i = too(2n + I) of an ideal gas in an isotropic harmonic potential. For the 
lowest radial modes with n — 0, the TF eigenfrequencies u>oVl differ from those for the harmonic oscillator ujqI except 
for the lowest dipole mode (1 = 1). 

The associated eigenfunctions are polynomials of the form n n i(r) = r l P n i(r 2 / R 2 ), where P n i(u) is a polynomial 
of order n (a type of Jacobi polynomial), for example, Pqi(u) = 1 and P\i(u) — 1 — (21 + 5)u/(2l + 3). The radial 
eigenfunctions are bounded at both the origin and the TF surface r = Rq. 



2. anisotropic axisymmetric trap 

The preceding treatment can be generalized to an axisymmetric but anisotropic trap with Vt r = \MiJ\(r\ + A 2 z 2 ), 
where A = u) z /u>± is the asymmetry parameter. The condensate density n is proportional to 1 — r^/R 2 ^ — z 2 /R 2 , 
where R± and R z are the radial and axial semiaxes of the ellipsoidal condensate with anisotropy ratio R z /Rj_ = 1/A. 



Equation (3.50) now has eigenfunctions n! proportional to e im ^ and eigenvalues that depend explicitly on the azimuthal 
quantum number m. They fall into two distinct classes 28 5^j5^| 



a. n' oc r ™ e im ^x a polynomial in (r^,z 2 ) (these modes are even in z); 

b. n' oc zrj"' e m ^x a polynomial in (r\,z 2 ) (these modes are odd in z). 

Numerous special cases are of interest, and it is simplest to start with the choice that the polynomial equals one. 

a. The solutions of type a take the form n' oc r l Yi m with m = ±Z, and the associated frequency is u 2 n =± l = lu>±, 
independent of the asymmetry parameter A. This mode is a generalization of that with frequency Wq; = lu 2 in a 
spherical trap. For 1 = 1, these two modes are the circularly polarized radial dipole oscillations with frequency u>± 
and m = ±1. 

b. The solutions of type b take the form n' oc r l Yi, n with m = +(l — 1), and the associated frequency is ^^=±(1-1) = 
(I — l)oi>± + u 2 . For I = 1, this mode is the axial dipole oscillation with m = and frequency u> z . 

The next simplest class takes a polynomial of the form A + Br\ + Cz 2 . For m = 0, those of type a have the two 
frequencies 



= 2 + f A 2 ± ^(2 - |A 2 ) 2 + 2A 2 . (3.53) 

The two eigenstates involve a coupled motion of a monopole mode with n = 1, I = and a quadrupole mode with 
n = 0, I — 2, m = 0. For a spherical trap (A = 1), the frequencies reproduce both the results given in eq. ( |3.52|) . For 
a disk-shape trap with A = \/8, this expression gives a good fit to the observations on 87 Rb |2(|. For a cigar-shape 

trap with A <C 1, the two frequencies reduce to io z and 2ui±; the first fits the observations on 23 Na [^7j within 

1%. More generally, theoretical eigenvalues and eigenfunctions are readily determined for larger values of \m\ and for 
the solutions of class b. 



E. Sum-rules for a non-uniform Bose condensate 

Stringari ]28fl has modified the sum rules for a uniform medium to obtain accurate estimates of the first few lowest 
frequencies of the collective modes in a trap for intermediate values of the interaction parameter Na/do, where neither 
the nearly ideal gas nor the TF approximation is accurate. Consider a general operator F that excites the system 
from its ground state |0) to an excited state \j ). Here, H\j) ~ Ej\j) and H\Q) — Eq\0), where H is the interacting 
many-particle Hamiltonian and Ej and Eq are the corresponding exact energies. 

In analogy with eq. ( 3.18| ), define a "spectral-density function" 



SM = £|0lF|0)| 2 5(c-^o), (3.54) 

where ojjo = (Ej — Eq)/Ti. In the following analysis, the matrix elements of F are assumed to have the following 
time-reversal behavior 

|0lF|0)| 2 = |0lFt|0)| 2 = |(0|F|j)| 2 , (3.55) 
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which restricts the allowed class of excitation operators. The various moments of the spectral-density function are 

p DO 

Mi= Ajjuj 1 S(u). (3.56) 
Jo 

For example, the first moment is 

Mi = ^2\(j\F\0)\ 2 uj j0 

= \ E (<°l^ f b' ) (j\F\0) + (0\F\j ) Uj0 (j\Ft\0] 



= ^E(< ^ t lj)0'l[^^l°>-< l[^^li>0'l^ t |0)). (3-57) 



The sum over the excited states can now be performed with the completeness relation [this is essentially the method 
used to prove the /-sum rule in eq. ( 3.25| )], leading to the expression 



F*,[H,F] |0), (3.58) 



expressed in terms of a ground-state expectation value of a double commutator. Similarly, 

1 



M 3 = n< 

2h 3 



[F*,H],[H,[H,F]] |0). (3.59) 



1. surface modes 

It is simplest to treat a spherical trap with a spherical condensate density. Consider the excitation operator 

N 

F = r\ Yi m (9i, fa) in first-quantized notation, or (3.60a) 

i=l 

F = J d 3 r ffi(r) r l Yi m (8, <fi)ip(r) in second-quantized notation. (3.60b) 

This operator has no radial node and thus predominantly excites modes with zero radial quantum number n = 
(called "surface modes" [p8[). A direct evaluation of the two moments for the ground-state condensate wave function 
(the solution of the GP equation) yields the ratio M3/M1, which is an upper bound for the squared frequency of the 
transition from the ground state to the first state excited by the operator F. In fact, this ratio should provide a good 
approximation to the exact lowest squared frequency since these surface modes are orthogonal to the higher radial 
modes. 

For the special case of the quadrupolc surface mode (n = and 1 — 2), this ratio gives the expression 

^ = 2 ^{ 1+ ^)' (3 ' 61) 

where the angular brackets are the ground-state expectation values of the kinetic and trap potential energies evaluated 
with the condensate wave function 'J. This expression holds for all values of the interaction parameter Na/dtj, 
providing a uniform interpolation between the ideal gas and the TF limit. 

a. For an ideal gas, the familiar properties of the harmonic oscillator show that the kinetic energy and trap potential 
energy are equal, reproducing the expected result that LU02 — 2loq for the excitation to the lowest quadrupole state of 
the oscillator. 

b. In contrast, the kinetic energy is negligible in the TF limit, reproducing the previous TF result from eq. ( 3.52| ) 
that W02 = V2i^o- 
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2. lowest compressional mode 



The operator 

N 

F = ^^r 2 , or, equivalently, (3.62a) 



i=l 



F = / d 3 rft(r)r 2 4>(r) (3.62b) 



excites the lowest monopole mod e with n — 1 and I — 0. Direct evaluation of the moments Mi and M3, along with 
the virial theorem from eq. ( [2.28 ) gives the result 

a. In the ideal-gas limit, the last ratio is 1, yielding the expected expression u>iq — 2ujq [from the general expression 
u> n i = Wo(2n + I) for an isotropic oscillator]. 

b. In th e TF limit, the last ratio is negligible, yielding coio = y5u>o, in agreement with the previous TF result from 
eq. ( p2l 



IV. VORTICES IN A DILUTE TRAPPED BOSE GAS 

The prediction and detection of quantized vortices in superfluid 4 He were early milestones in understanding the role 
of the condensate and the associated macroscopic wave function (see, for example, refs. [|35|j56| ). Thus it is natural to 
consider similar vortex structures in the dilute trapped Bose gases, but the analysis involves several new and subtle 
features. 



A. Review of classical vortices 



Consider an incompressible nonviscous fluid with uniform number density n. A long straight vortex line at the 
origin has the velocity field 

v(r) = (4.1) 

where <j> is a unit tangential vector in cylindrical polar coordinates. The flow velocity consists of circular streamlines, 
and its magnitude diverges near the symmetry axis (as rj_ — > 0). It is obvious from the flow pattern that V • v = 0, 
so that the flow is transverse (or solenoidal). In addition, it is not difficult to verify that V x v = except on the 
axis, and use of Stokes's theorem shows that 

Vxv = di (2) (ri). (4.2) 

Thus the flow is longitudinal (or irrotational) nearly everywhere, but it is singular on the axis (analogous to the 
magnetic field of a long straight current-carrying filament) . Finally, the line integral of v is known as the "circulation;" 
for any contour C that encircles the vortex once, the circulation 

dl-v = K, (4.3) 

c 



is independent of C, with the constant k in eq. ( J4.1| ) fixing the strength of the vortex. 

The energy is purely kinetic in an incompressible fluid. In the present case of a long vortex filament, it has the 
value (per unit length) 



E v = \Mn J d 2 r± _ v 2 



MfiK 2 f dr± MhK 2 ( R \ 

In — . (4.4) 



Here, the integral must be cut off at an upper limit i?o, which can be interpreted as either the size of the container 
or the intervortex separation, and a lower limit £, which can be interpreted as the size of the vortex core. 
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1. quantized circulation 



The idea of a quantized vortex in superfluid He can be motivated by recalling the Bohr-Sommerfeld quantization 
condition § d \ ■ p = h for the periodic motion of a particle in phase space. Since p = Mv for a particle in the vortex, 
this result immediately suggests that the circulation in superfluid 4 He should be quantized in units of h/M. Note that 
k has the dimensions of length 2 /time (like a diffusion constant or a thermal conductivity or a kinematic viscosity). 
Onsager |57| first publicly proposed the concept of quantized circulation in 1949 (sec. 2.3 of rcf. contains a brief 
history of Onsager's unpublished work). London |58| documented one of Onsager's earlier remarks and made a similar 
proposal (59] concerning quantized magnetic flux in superconductors in 1950. Independently, Feynman developed a 
more complete description of quantized vortex lines in superfluid 4 He f)0"| , based on the identification of the velocity 
with the gradient of the phase of the appropriate many-body wave function [compare eq. ( |3.36 ) for the condensate 
velocity] . 

A classical vortex in an incompressible fluid requires a specific model for the core radius £. In fact, all real fluids 
have a non-zero compressibility, implying a finite speed of sound s. If the hydrodynamic flow speed |w| is small 
compared to s, then the fluid is essentially incompressible. This picture necessarily fails near the vortex core, where 



eq. (4.1) shows the the flow speed diverges; hence a singularity develops when the magnitude of the circulating velocity 



becomes comparable with the speed of sound. Setting v(£) = k/27t£ « s gives the resulting core radius 

«"S5- (4 - 5) 

If the vortex is quantized with circulation h/M, this result implies that the vortex core has a radius 

(4-6) 

Ms 



As seen in eq. (2.43), the same qualitative expression holds for the healing length £ = h/(y/2Ms) in a dilute trapped 



Bose gas in the TF limit, suggesting that £ also characterizes the vortex core radius in such a condensate. 



2. classical vortex dynamics 



The dynamics of a vortex in an incompressible fluid follows from the celebrated "circulation" theorem of classical 
hydrodynamics (see, for example, sec. 48 of ref. Q]); it states that the circulation around any contour that moves with 
the fluid is a strict constant of the motion. If the contour initially contains a vortex with circulation «, the vortex can 
never escape the moving contour. In the particular case of a single long straight vortex line in otherwise stationary 
fluid, this result implies that the vortex itself remains at rest. Similarly, in the presence of a uniform external flow 
v GX , the vortex line moves rigidly with the same velocity (this conclusion also follows with a Galilean transformation 
to a moving coordinate system). 

A more interesting case is two long straight parallel vortices a distance d apart. Since neither vortex moves under 
its own influence, the motion arises solely from its neighbor. Assume first that they have the same sense of circulation, 
as shown in fig. 5a. In this case, it is easy to see from the circulating velocity fields that they each execute a circular 
trajectory at fixed separation d. Since the interaction energy depends only on d (and the individual circulation k), 
the constancy of the separation d follows from the conservation of energy in a nonviscous fluid. In the presence of 
weak dissipation, the system will act to lower the energy; here, the two vortices will slowly spiral apart because the 
total energy of two identical vortex lines is lowest at infinite separation (for d — > 0, the vortices fuse into a single 
vortex with circulation 2k and four times the energy of a single vortex line, whereas, for d — ► oo, the energy reduces 
to two times the energy of a single vortex line). 

If the parallel vortices have opposite sense of circulation, they are known as a "vortex pair" (shown in fig. 5b). 
The mutual influence of the circulating velocity fields means that they move perpendicular to the line joining them 
in the same sense as the local fluid between them, with fixed separation d (again because of conservation of energy) . 
This configuration is the two-dimensional analog of a vortex ring, which moves perpendicular to its plane maintaining 
a fixed radius. In the presence of weak dissipation, the two vortices slowly drift together and annihilate when the 
separation d is comparable with the vortex core size £. 



3. effect of rigid boundaries in classical hydrodynamics 

As a simple example, consider a long circular cylinder with rigid walls of radius Rq. If a vortex line is located on the 
symmetry axis, its circulating flow velocity is everywhere parallel to the boundary, so that the flow field is unperturbed 



31 



by the presence of the cylinder. The situation is more complicated, however, if the vortex is shifted rigidly a distance 
ro < Ro off-axis. In this case, the flow of the vortex by itself no longer matches the boundary condition that the normal 
component of the flow velocity must vanish. This boundary-value problem has an elementary solution consisting of 
a single opposite image located along the same ray at a distance Rq/to outside the cylinder (see fig. 6). Comparison 
with fig. 5b shows that the vortex and its image constitute a vortex pair; in the absence of dissipation, the original 
vortex will therefore execute a circular orbit at fixed ro under the influence of its image (which moves synchronously 
with the vortex). If the system is weakly dissipative, the vortex will slowly spiral outward and eventually annihilate 
with the image vortex when tq « Rq — £. 



4- effect of rotation 



One of the original motivations for the introduction of quantized vortices was the observed behavior of a rotating 
superfluid. For example, a bucket of 4 He rotating at an angular velocity fl was known [5^ ] to retain its classical 
parabolic meniscus z = ^il 2 r 2 /g far below the superfluid transition temperature T\. This observation appeared to 
contradict Landau's assertion that the superfluid velocity is irrotational with V x v s = (implying that it could 
not rotate) and that p s / p ~ 1 (so that the normal fluid played a negligible role). To explain this apparent paradox, 
Feynman J60| ] relaxed the irrotational condition of strictly zero vorticity, allowing the vorticity to be singular at the 
cores of the quantized vortices [compare eq. (4-2)]; in particular, he suggested that a rotating superfluid contains an 
array of quantized vortex lines with areal density n v — 2Q/k = MQ/(ttJi) and mean vorticity n v n — 20, thereby 
mimicking the solid-body rotation required to explain the curved meniscus (see sec. 2.4 of ref. p6fl ). 

^From this perspective, it becomes interesting to consider the effect of rotating a physical container at an angular 
velocity ft — 0,2. In the laboratory (non-rotating) frame, the moving walls represent time-dependent potentials that 
do work on the system, precluding a description in terms of thermodynamic equilibrium. In the rotating frame, 
however, the walls are stationary, so that the hamiltonian becomes time-independent, and the usual thermal Gibbs 
distribution remains valid when expressed in terms of the energy levels appropriate for the rotating frame. In the 
rotating frame, the new hamiltonian is known to be H — ■ L = H — flL z (see, for example, sec. 5 of ref. jll) and 
sec. 34 of ref. [§). 

At zero temperature, a system rotating at an angular velocity f2 will adjust any free parameters to minimize the 
quantity F = (H — ilL z ) = E — flL z , which can be considered a "free energy" that depends on the parameter SI. As a 
simple example, consider a long cylinder of radius Rq- If the fluid does not contain a vortex, the velocity is everywhere 
zero, so that the energy Eq and angular momentum Lq per unit length both vanish, with Fq = Eq — CILq = 0. When 
a vortex is placed at the center of the cylinder, however, the energy per unit length is 



Mm 2 (Rq 
E v = — In — 



(4.7) 



from eq. (4.4), and an elementary calculation shows that the corresponding angular momentum per unit length is 
L v = ^MukRq. Thus the free energy per unit length of the fluid containing one vortex is given by 



F„ — E„ — flL„ — 



-in 



In I -j- ) - \nMnnRl. 



(4.8) 



This quantity is positive for sufficiently small f2, but it decreases linearly with increasing angular velocity and becomes 
negative (and hence less than Fq for the no-vortex state) at a critical angular velocity O c i = (k/2ttRq) ln(i?o/£) f° r 
the creation of a single vortex line in a long cylinder of radius Rq. Use of the quantized circulation k = h/M yields 
the expression (see, for example, sec. 2.4 of ref. |56|) 



AIRq 



111 



Rq 



for a quantized vortex in uniform fluid. 



(4.9) 



This elementary analysis can be generalized to the case of a long straight vortex a distance rp from the center of 



the cylinder. The free energy per unit length remains zero if there is no vortex, and a detailed calculation 
gives the free energy per unit length of a vortex displaced a fractional distance xq — tq/Rq from the center 
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(4.10) 
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where fio = k /(2tt^o) = ^/{MRq) is a characteristic angular velocity. Figure 7 shows this quantity as a function of 
xq for various fixed values of f2. If = 0, the free energy decreases monotonically with increasing xq, indicating that 
the vortex is unstable in the presence of a weak dissipation. This behavior persists to the value fl — f^o, when the 
curvature at the origin xq = vanishes. For f^o < ^ < ^ci, the free energy develops a local minimum at xq = 0, 
although F v (0) > F v (l), with an intermediate barrier that hinders the vortex from reaching the outer wall at xq = 1. 
This behavior suggests that f2 m = Qq represents the onset of metastability in the presence of weak dissipation, which 
persists up to the critical angular velocity fid, when the free energy F v (0) for a vortex at the origin becomes equal 
to that F v (l) for a vortex at the wall. For Q, c i < f2, the vortex becomes stable relative to the no- vortex state. 
Experiments on rotating superfluid 4 He p2| have verified these predictions in considerable detail. 



B. Vortices in a dilute uniform Bose gas 

The original application of the GP equation pO|pl| was the description of a long straight vortex in an otherwise 



uniform dilute Bose gas. Recall eq. (3.36) for the condensate velocity vo = (?i/M)V5, where S is the phase of 
the condensate wave function, and consider the circulation T = <f c dl-vo around an arbitrary closed contour C. A 
combination of these expressions gives 

r = Ai d i.vs = A A 5L l (4.ii) 



M J c M lc " 

where AS\c is the change in the phase on once going around C. If the condensate wave function is single valued, this 
change must be an integer times 2tt, implying that the circulation in any condensate is quantized in units of h/M . 
In general, the GP equation for a bulk dilute Bose gas is 

(T + .g|*| 2 - A i)* = 0, (4.12) 

since Vjr is absent. Assume that the condensate contains a singly quantized vortex line at r± — 0, oriented along z, 
with bulk condensate density no far from the vortex (r± — > oo). The condensate wave function has the form [ p0]]2l| ] 

* = Vn?e*/(rx), (4.13) 

where / is real and approaches 1 as rj_ — ► oo. Here, the phase is given by S — cj>, so that the condensate velocity 
vo = (h/M)V(f> = (Ti/Mr±)(f> is precisely that for a vortex with circulation h/M. The condensate density |^| 2 = 
n o( r ±) — n of{ r ±) 2 is cylindrically symmetric. The radial amplitude obeys the nonlinear GP equation 

h 2 1 d d h 2 
■ r 



2Mr ± drj_ dr ± 2Mr 2 , 



gnof - M / = 0, (4.14) 



where the asymptotic behavior of / fixes the chemical potential a = gn at the value for the uniform system. 

In a dilute uniform Bose gas, the healing length is defined by £ 2 = Ti 2 / (2Mgno) = l/(87rano), balancing the kinetic 
energy and the interaction energy. With the dimensionless variable x — r±/£, the radial GP equation becomes 

-^ x ^ + -2+f 2 - l )f = ^ ( 4 - 15 ) 

x dx dx x z J 

where the term 1/x 2 is the centrifugal barrier associated with the circulating velocity field. For small x <C 1, this 
barrier forces the amplitude to vanish linearly inside the vortex core x ^ 1, and a detailed analysis shows that 

tf r \~ f 0-583 x, forx<l, , 

W~ \i-(2x 2 )-\ forx»l. (4 - lbj 

Figure 8 shows the form of this function, indicating that this quantum vortex autom atica lly forms its own core with 
radius w £ = h/(\/2Ms), where s is the Bogoliubov speed of sound [compare eq. fl4.6D l- The condensate density 
^o(^_l) varies like r\ near the origin, and the corresponding condensate current jo(?\L) = "■o(' , ±)vo(? , j^) reaches a 
maximum near the core and then vanishes linearly as r± — > 0, even though the velocity itself diverges. Numerical 
analysis yields the GP vortex energy per unit length |^3| 

E v « _ In (^1.46— j « ^— In (^1.46- j , (4.17) 
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where Rq is an outer cutoff; this result has precisely the same form as the classical expression in eq. (4.4). 

The GP description of a vortex in a bulk dilute Bose gas can be generalized to study the dynamics of a set of widely 
separated parallel vortex lines |64|] located at {rj_j}, assuming that |rj_j — r±j\ ^> £ for all pairs. An approximate 
GP condensate wave function can be constructed as a product of the individual radial functions /(|rj_ — v±j\), each 
of which approaches 1 far from the vortex cores; the phase is given as the sum of the azimuthal angles Sj measured 
relative to r±j as origin, where tan Sj = (y — yj) / \x — Xj), so that 

*(r ± ,t) = V^~oe-^ t/h l[[e^f(\r ± - r ±j \)] . (4.18) 

j 

This state changes with time both because of the overall chemical potential and because each vortex follows its own 



dynamical trajectory. Substitution of this condensate wave function into the time-dependent GP eq. (2.24) determines 
the subsequent motion of the individual vortices; they obey the classical hydrodynamic prescription that each vortex 
moves with the local velocity induced at its position by all the other ones. This result should not be very surprising, 
because the GP equation itself can be recast in a hydrodynamic form. 



C. A vortex in a dilute trapped Bose condensate 



Assume a singly quantized vortex in an axisymmetric trap, where the condensate wave function takes the form 

9(T) = y/Noe**f{r ± ,z). (4.19) 

Like the case of a vortex in a uniform condensate, the condensate velocity field here is just Vo = (h/Mr±)<f>, represent- 
ing a vortex with circulation k — h/M . The amplitude function / satisfies the GP equation and can be determined 
numerically [p2f . This function vanishes linearly as r± — > 0, giving a node on the symmetry axis. Consequently, the 
condensate density for a vortex in a trapped condensate is toroidal (in contrast to the non-vortex condensate, where 
the density decreases monotonically away from the center). 



1. critical angular velocity Q, c i 



Given the condensate wave function for the vortex state and that for no vortex, it is straightforward to evaluate 
the additional energy AE associated with the formation of the vortex, along with the angular momentum L z . If the 
system rotates with angular velocity Q, a vortex line on the symmetry axis becomes stable when AF = AE — flL z 
becomes negative, so that the critical angular velocity is fi c i = AE/L Z . 

This quantity can be evaluated analytically in two limits. If the system is nearly ideal, with Na/do <C 1, then the 
condensate wave function for the vortex is one with all No particles in the single-particle state 

Xio(t±)M*) « rj_e* exp (-^ - , (4.20) 

where x«+n_ is a two-dimensional oscillator state with n+ and n_ right- and left-circular quanta and i\> n is the usual 
one-dimensional oscillator wave function J65| . The increased energy is essentially N times the increased single-particle 
energy hoj± associated with the transition from the two-dimensional ground state %oo to the vortex state xw, an d 
the associated angular momentum is just NqH because every particle has unit angular momentum. Thus fl c i — > lu± 
for an ideal gas in an axisymmetric trap. In this ideal-gas limit, however, the transition is hugely degenerate, for the 
same f2 c i also characterizes a condensate with multiple quanta of circulation. With the GP equation, it is not difficult 
to evaluate the first correction to this value for a singly quantized vortex in a nearly ideal Bose gas 

J= ^ + -.., (4.21) 

showing that the critical angular velocity decreases linearly as the interaction parameter Noa/d z = \/XNoa/dx 
initially increases from (here, A = lo z /lo±_ is the asymmetry parameter for the trap). 

In the opposite l imi t of a large interaction parameter Na/do 3> 1, initial estimates p5j relied on classical hydrody- 



namics, using eq. (4.9) to write 
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fLi ft , , (R 



J rln ^ =^m(^], (4.22) 



since = h/Mu±. Note that this ratio is small in the TF limit where d\jR\ ~ £,/R± <C 1; the explicit TF relation 

2/5 



dj _ / d ± V 



(4.23) 



shows precisely how this small parameter scales with TV and A. 

An improved TF estimate follows by writing the condensate wave function as = e 1 ^ I \j/ 1 . and ignoring the gradients 
of \&\. The new feature is that the TF density profile now has a centrifugal barrier [pqpTfl 



h 2 \ ( h 2 



gni « U _ |Mo;iri - \Mu? z z 2 - j 6 (mj - iM^iri - |M W ^ - , (4.24) 

where ni and /xi are the density and chemical potential f or th e state with one quantum of circulation. This density 



can be compared with that for the no-vortex state in eq. (2.35) 



gn « (/i - \Muj\r\ - \Mu 2 z z 2 ) 9(/i - \MuJ\r\ - \Muj\z 2 ) . (4.25) 

The fractional change in the chemical potential caused by the vortex (hi — ^io)/Mo can be shown f68|| to be small, of 
order (d^_/R^_) l^R^/d 2 ^. As a result, the TF density profile for the condensate with a singly quantized vortex has 
the simple form 



Rl Rl r 2 J U { L R\ Rl J 



m(r ± , z) « no(0) 6 1-^-^-^ , (4.26) 



where no(0) is the central density for the no- vortex condensate. The important new qualitative feature of a vortex 
in the TF limit is the appearance of a small flared hole of radius ~ £ [ ^2p^j67|] , but the remainder of the condensate 
density is essentially unchanged; in practice, it is usually sufficient to retain the no-vortex density and simply cut off 
any divergent radial integrals at the core size £. 

A detailed solution of the GP equation for an axisymmetric trap in the TF limit gives the more careful 

estimate 

5ft, (R±\ . , , fi c i 5 d 2 /R ± \ /Ann ^ 

^ lK 2MRl ln \—) OT > e ^ vale ^ — "2i?i ln vTj' (4 ' 27) 

which holds with logarithmic accuracy; it is larger than the previous naive estimate in eq. (4.22) by a factor I 



that arises from the reduced angular momentum for the non-uniform density relative to that for a uniform density. 
Reference |34| shows that this asymptotic estimate agrees well with the numerical work of ref. |^2| for the particular 
asymmetry A = \/8, appropriate for the original experiments juj on 87 Rb. 

The TF description has also been extended p9[ | to the case of a rotating axisymmetric trap containing a straight 
vortex that is displaced laterally a distance r$ from the center. The dependence of the free energy on the dimensionless 
parameter xq = r$/R± for fixed 17 is qualitatively similar to that in fig. 7, but the details differ consi derab ly. In 



particular, the vortex at the center is unstable for f2 < f2 m = |^ c i, where £7 cl is given above in eq. (4.27), and 
metastable for fl m < O < f2 c i. Thus the presence of the trap reduces (but does not eliminate) the regime of 
metastability. A similar result also holds for a simpler classical model where the uniform density n is replaced by a 
parabolic radial density profile {nj. In both these cases, the reduced outer density in a trap lowers the free energy 
F v {xq) for xo ^ \ compared to that for uniform density, raising the threshold for metastability relative to that in 
fig- 7. 



2. possible scenario for creation of a vortex 

A basic question is: how can the condensate be made to rotate? For 4 He in a rotating circular container, the 
microscopically rough walls can spin up the normal (viscous) fluid when the temperature is above the superfluid 
transition. If the fluid is then cooled below the transition temperature, the superfluid is created in a state of rotation 
for sufficiently large rotation speeds. This method fails for a condensate in an axisymmetric trap, because the trap 
is simply a smooth potential. Consequently, only a non-axisymmetric trap can spin up the condensate, for the 
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asymmetry in the trap potential exerts a torque on the condensate, similar to the effect of moving walls in classical 
hydrodynamics. 

It is helpful first to consider the classical example of a uniform fluid in a long rotating elliptical cylinder with semi- 
major and semi-minor axes a and b. If a > 6, the rotating walls push the fluid, creating a classical velocity potential 
<I> c i = flxy(a 2 — b 2 )/(a 2 + b 2 ) [7l|, even in the absence of a vortex. The resulting angular momentum L z q = IoCl 
and kinetic energy Eq = ^Iq£1 2 have the expected dependence on the rotation speed f2, with an effective moment of 
inertia per unit length Iq that is less than the classical solid-body value J s b by the ratio 

(4.28) 



/ sb \a 2 + b 2 

This quantity vanishes for a circular cylinder (a = b), but it can approach 1 from below for extreme asymmetry 
(i < a). Unlike the case of a circular cylinder, the free energy Fq = Eq — HL z o = —^Iott 2 for a vortex- free state is 
now negative. Note that the classical velocity v c i = V4> c i is everywhere irrotational with zero vorticity V x v c i = 0. 
Nevertheless, the moment of inertia Jo is non-zero; indeed, for b -C a, it can be comparable with that for solid-body 
rotation v s b = ^2 x r, whose vorticity is uniform with V x v s b = 2J7. 

The critical angular velocity f2 c i for creating a vortex at the center dep ends on the asymmetry ratio b/a |73], and 
experiments confirm the theoretical predictions in considerable detail |73j. In the extreme asymmetric limit b -C a, a 
detailed calculation yields 




2Mb 2 V£ 

which holds with logarithmic accuracy (here, the circulation is taken as h/M). Qualitatively, this value can be 
interpreted as that for a circular cylinder with the radius equal to the largest inscribed circle that fits inside the 
elliptical cross section. 

An approximate TF solution of the GP equation in a rotating triaxial (totally asymmetric) trap with to x =/= lu v lu z 
bears out these qualitative conclusions ]69[ |. As usual in the TF limit, the chemical potential fi determines the 
condensate radii i?„ = 2[ijMuj 2 a . Even in the absence of a vortex, the rotating walls push the condensate, inducing a 
phase for the condensate wave function that is linear in the angular velocity 1~2 = ilz 

This expression has the correct angular symmetry oc xy = ^r 2 sin 20 and vanishes for an axisymmetric trap; it is the 
quantum analog of the classical velocity potential for a rotating elliptic cylinder |7lj| . 

For a non-axisymmetric trap rotating about the z axis, it is convenient to define the mean radius 

a 2R l R l ,,211 , 

^ = if+if' o r i equivalently, ^ = ^ + ^ (4.31) 

so that 1Z± — * = R x = R y for an axisymmetric trap. A calculation of the energy and angular momentum of a 
vortex on the axis of the trap yields the critical angular velocity p9[ | 

for the creation of a vortex in an arbitrary non-axisymmetric trap in the TF limit. For an axisymmetric trap 



(1Z± — > R±), the resulting expression for f2 c i reproduces the previous result in eq. (4.27); in contrast, it becomes 

R 
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in the extreme anisotropic limit R y -C R x , very similar to the classical expression in eq. ( 1.29| ) 



One possible experimental scenario for creating a vortex in a trap follows from the analogy with superfluid 4 He 
]62[ . Cool a rotating trap with significant anisotropy in the xy plane from the normal state through T c to some low 
temperature T <g; T c . If fl exceeds the critical angular velocity O c i, the resulting Bose condensate should have a 
vortex. Note the crucial order of operations, with the condensate created in a rotating state. Although reversing 
the order (cooling and then rotating) does often create a vortex in superfluid 4 He, such a procedure is hysteretic and 
displays considerable metastability J74|;[75|]. Type-II superconductors exhibit similar metastability when a magnetic 
field is applied after the sample has been cooled in zero field. 
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3. possible detection of a vortex line 



In principle, it may be possible to observe directly the reduced density in the vortex core (somewhat like the optical 
detection of localized trapped charge in the core of vortex lines in superfluid 4 He ) , but the small healing length 

£ in the TF limit may well suppress the effect. An alternative and promising approach considers the effect of a vortex 
on the collective modes of the condensate. The basic physical idea is that a uniform rotation splits the time-reversal 
degeneracy of the two modes with ±|m| that are originally degenerate (like the Zeeman effect of a uniform magnetic 
field on the states with different azimuthal components of magnetic moment). 

Two distinct approaches have been developed (they both yield the same results in the TF limit). One uses sum 
rules j76| to generalize the earlier study of the low- lying collective modes for a vortex- free condensate |2q] , discussed 
in the previous section. The other uses the hydrodynamic equations in the TF limit [jS8|| , where the vortex provides 
a weak perturbation (ref. |56| also obtained some, but not all, of the same results). 

To estimate the magnitude of the effect, recall that a static velocity field Vo alters the usual time derivative d/dt 
to the hydrodynamic derivative d/dt + vo -V. Thus the original frequency u) of a wave is shifted to ~ to ± vq/Rq, 
where Rq is the typical dimension of the condensate. For a vortex, the circulating velocity is of order vq ~ %/ {MRq), 
and the fractional change in the frequency is 5uj/lu ~ %/(MRqU>±) ~ d\/R\ <C I in the TF limit. Alternatively, the 
fractional change in the frequency is comparable to the ratio vq/s, where s is the speed of sound. In a dilute trapped 
Bose gas, this speed is of order s ~ Ti/(M£), which becomes s <~ T%Ro/(Md l ) in the TF limit and reproduces the 
previous estimate for 5u>/u>. Similar ideas (see, for example, pp. 144-145 of ref. pqj ) have served to study the effect 
of a quantized vortex line in a type-II superconductor on the BCS quasiparticles. 

Consider a large (TF) vortex-free condensate in an axisymmetric trap, and focus on a particular collective mode 
with degenerate frequency u? m ,. In the presence of a vortex, the fractional splitting of the modes with ±|m| can be 
characterized by the ratio 



H - ^-H _ . d\ _ A / dx \ 

fi — L\\m.\ T^T ~ TZTTZ — T 



W ... 



2/5 



where Ai m i is a dimensionless number that depends on the specific mode in question. For \m\ > 1, it is proportional 
to the matrix element of 1/tj_, evaluated with the normalized unperturbed (vortex- free) eigenf unctions, but the case 
of \m\ = 1 requires a more careful analysis j pSfl . 

a. For the lowest quadrupole mode with |m| = 2, the unperturbed frequency is ui® = V^^x, and the dimensionless 
factor is A 2 = 7/y/2. 

b. For the lowest quadrupole mode with |m| = 1, the unperturbed frequency is u>® = vl + A 2 lux = \J &x + u>\ , 
and the dimensionless factor is Ai = 7A 2 /(1 + A 2 ) 3 / 2 . 

Typically, the predicted fractional splitting is of order 10% for the low-lying modes. Current measurements of 
frequencies for non- rotating condensates are accurate to % 1%, so that this vortex-induced splitting should be readily 
detectable. 



4- stability of a vortex in a trapped Bose condensate 

If the condensate rotates with an angular velocity il > il c i , a singly quantized vortex is believed to be the stable 
equilibrium configuration. If, however, the angular velocity is reduced so that £1 < f2 c i, does the vortex remain stable 
(especially if f2/f2 c i <C 1)7 In the present case of a trapped condensate, the question can become even more intricate, 
depending on whether the trap remains non-axisymmetric (so that the reduced rotation exerts a braking torque on 
the condensate) or is first adiabatically transformed to axisymmetric (in which case, the reduced angular velocity of 
the trap potential cannot affect the rotation of the condensate). 

Rokhsar |67| has given a qualitative argument that a vortex in a non-rotating axisymmetric nearly ideal condensate 
should be unstable. In this limit, the single-particle vortex wave function has the form yj v w e 4 ^|^„|, with the radius 
of the vortex core comparable to the oscillator length and condensate radius do (note that the coherence length £ 
here is much larger than do and therefore does not characterize the vortex core radius in the nearly ideal limit). This 
large depleted region of reduced repulsive Hartree potential energy favors the formation of a localized bound core 
state, with a real single-particle wave function ?/V The actual condensate wave function becomes a linear combination 

« y/N v ipv + yNbipb, where N v is the number of particles in the vortex condensate and Nb = N — N v is that in 
the bound core state. For small Nf,, this complex wave function w \/N v e' l '''\ip v \ + ^/Nb\ipb\ has a node away from 
the axis of symmetry, so that the vortex core shifts radially outward from the origin. This motion can occur only if 
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there is some mechanism for the condensate to lose energy and angular momentum; it is a quantum analog of the 
hydrodynamic instability below the threshold for metastability, discussed in connection with fig. 7. As in that case, 
the time for a vortex to move outward depends on the details of the dissipation. In addition, the vortex core itself 
becomes small in the TF limit and cannot support a bound state for Na/do > 1 f67j| , so that this argument is not 
applicable in the TF limit. 

As an alternative physical picture in the non- interacting limit, note that the vortex has all the particles in the 



single-particle state xio(r±)ipo{z) given in eq. (4.2C), with excitation energy Hlo± and angular momentum h. Clearly, 
the system can lower its energy by transferring one particle from this (excited) condensate to the original ground 
state xoo( r ±)'4'o(z); the same process can be repeated many times, which is the physical source of the instability. 

This argument can be sharpened somewhat by considering the Bogoliubov equations for a singly quantized vortex in 
a nearly ideal condensate with Na/d z <C 1 |7Cj ]. Among the many normal modes with positive norm J d 3 r(\u\ 2 — \v\ 2 ) = 
1, only one is anomalous and has a negative frequency 

LU„ 1 And 

— ,— , +■■■■ 4.35 

w± V8n d z 

This mode consists of a coherent linear combination of both u and v amplitudes with an angular dependence m a = — 1 
relative to the vortex. When re-expressed in terms of the particle density, it represents a displacement of the vortex 
core relative to the center of mass of the condensate. 

More generally, Dodd et al. ^fi\ have studied the normal modes of a singly quantized vortex numerically for small 
and intermediate values of the dimensionless coupling parameter Na/do (note that ref. |77j uses m to denote the 
quantum number of the static vortex condensate, whereas m here denotes the angular dependence oc e zm ^ of a density 
perturbation |6q| , in direct analogy with the same m dependence of a spherical harmonic Yj TO ). For the lowest-lying 
dipole modes with |m| = 1 relative to the vortex, Dodd et al. |77j find two dipole-sloshing modes with frequency luj_ in 
which the condensate oscillates rigidly with right- and left-circular helicity, independent of the interaction parameter 
(as expected). In addition, they exhibit an additional (anomalous) dipole mode with m = — 1 that has the frequency 



luj_ for the ideal gas and decreases toward zero with increasing coupling strength. Rokhsar |67| argues that this mode 
actually describes an antiparticle; it is the opposite-frequency and negative-norm partner of a physical mode with 
negative frequency and positive norm. With this reinterpretation, the frequency of the anomalous mode reduces to 
the expected value —uj± for an ideal gas and grows toward zero from below with increasing Na/d z . In the context of 
the preceding discussion, this anomalous mode is precisely the unstable (negative-frequency) mode found analytically 



in eq. (4.35) from the Bogoliubov equations in the weak-coupling limit. At present, little is known of its behavior in 
the TF limit, when Na/d z ^> 1, although the trend of the numerical data (77) suggests that the frequency remains 
negative. 

As noted several times, an applied rotation can stabilize a vortex. Specifically, the frequency tjj of a given normal 
mode in a non-rotating condensate is shifted to ujj — flrrij when the system rotates with angular velocity Q. In 
particular, the anomalous mode with negative frequency u a has m a = — 1; it therefore has an apparent frequency 
uj a + Q when observed in the rotating frame, indicating that a positive rotation raises the (negative) frequency, tending 



to stabilize the mode. Recall that eq. (4.21) gives the thermodynamic critical angular velocity f2 c i for the creation of a 



singly quantized vortex in the near-ideal limit, and comparison with the independently computed negative frequency 



uj a in eq. (4.35) shows that uj a + O c i vanishes identically for small Noa/d z . Hence, a rotation at tt c i just suffices 
to stabilize this unstable anomalous mode in the weak-coupling limit. The corresponding situation in the TF limit 
remains unknown, where the metastable frequency Q m < f2 c i discussed in connection with fig. 7 may well be relevant. 
In principle, various other methods can stabilize a vortex in a dilute trapped Bose condensate; in the present 



context, the most relevant proposal is to force the non-rotating condensate to assume a toroidal form [j78|j79| (for 
example, piercing it with a blue-detuned narrowly focused laser beam that acts as a repulsive localized force). In this 
case, the externally imposed hole in the stationary condensate renders it multiply connected, even in the absence of a 
vortex. Strictly speaking, the excited states of interest here are those with irrotational flow and quantized circulation 
around the circumference of the torus, but they are frequently considered to represent the corresponding quantized 
vortex with its core "pinned" in the central void. Such states have been studied in detail for superfluid He in a 
rotating annulus (see, for example, sees. 2.6 and 5.3 of ref. |56|]), where a row of physical vortices eventually appears 
in the gap between the walls at sufficiently high angular velocity. 
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FIG. 1. Dime nsion less excitation-energy spectrum E/ng as a function of the dimensionless wavenumber (a) Bogoliubov 



spectrum in eq. (1.42); (b) free-particle spectrum. 



FIG. 2. Bogoliubov coherence factors (a) u and (b) v from eq. (1.37) as functions of the dimensionless wavenumber fc£ 



FIG. 3. Dimensionless squared Bogoliubov energy E 2 /n 2 g 2 in eq. ( 2.44 ) for a uniform system subject to an attractive 
interaction with s-wave scattering length —\a\, as a function of the dimensionless wavenumber k£, where 87rn|a|£ 2 = 1. In these 
units, the unstable range is < fc£ < 



FIG. 4. Dimensionless variational energy 2E/(Nhu>o) in eq. ( 2.47) for a trapped condensate, as a function of the dimensionless 



parameter /3 that characterizes the actual radius /3do for various values of the interaction parameter: (a) Na/do = 0.5 (locally 
stable), (b) Na/do = 0.67 (onset of local instability), and (c) Na/do — 0.84 (locally unstable). 
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FIG. 5. Induced local velocity for two parallel vortex lines a distance d apart in unbounded incompressible dissipationless 
fluid: (a) Two parallel vortex lines with same sense of circulation; induced motion as shown yields circular orbits at fixed 
separation d; (b) two antiparallel vortex lines with opposite sense of circulation; induced motion as shown yields uniform 
motion perpendicular to line joining centers in direction of fluid at center, maintaining fixed separation d. 

FIG. 6. Cylinder of radius Ro containing a long straight vortex line displaced from the center by a distance ro < Ro- One 
image vortex with opposite circulation on same ray at the exterior point with radius Ro/ro suffices to satisfy the boundary 
condition that the normal component of velocity vanish at the boundary. 



FIG. 7. Dimensionless free energy per unit length MF v /(irnh 2 ) from eq. ( 4. 1C ) for a vortex in a cylinder of radius Ro; the 
vortex is displaced radially a fractional distance a;o = ro/Ro from the axis of symmetry, and the different curves correspond to 
various fixed values of the angular velocity O: (a) SI = (unstable); (b) 57 = Q m = f?o = ?i/(MRq) (marginally metastable at 
origin); (c) Q = 0.5f2 c i (metastable at origin); (d) = Q c i (stable at origin), where Q cl = f7 ln(i? /£) = (K/MRl) m(iio/£) 
is evaluated for Ro/£ ~ 100. 



FIG. 8. Radial wave function f(r±/£) obtained by numerical solution of eq. (4.15) 
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